Spatially Resolved Near-Infrared Spectroscopy of Seyfert 2 
Galaxies Mk 1066, NGC 2110, NGC 4388, and Mk 3 1 
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ABSTRACT 

We present near-infrared spectra with resolutions of A/AA ~ 1200 in the emission lines Pa/3, 
[Fe II] (A=1.2567/im), Br7, and H2 v=l— 0S(1) of the nuclei and circumnuclear regions of the 
four Seyfert 2 galaxies Mk 1066, NGC 2110, NGC 4388, and Mk 3. All of these galaxies show 
strong near-infrared line emission that is detected at radii several times the spatial resolution, 
corresponding to projected physical scales of 0.07 to 0.7 kpc. Velocity gradients are detected in 
these nuclei, as are spatial variations in line profiles and flux ratios. We compare the spatial and 
velocity distribution of the line emission to previously observed optical line and radio emission. 
The evidence indicates that the [Fe II] emission is associated with the Seyfert activity in the 
galaxies. Our data are consistent with X-ray heating being responsible for most of the [Fe II] 
emission, although differences in [Fe II] and Pa/3 line profiles associated with radio emission 
suggests that the [Fe II] emission is enhanced by fast shocks associated with radio outflows. The 
H2 emission is not as strongly associated with outflows or ionization cones as is the emission in 
other lines, but rather appears to be primarily associated with the disk of the galaxy. 



1. Introduction 

The near-infrared wavelength range offers valu- 
able probes of the dynamics, composition, and 
excitation mechanisms present in the circum- 
nuclear gas of active galactic nuclei (AGN). 
Among the near-infrared emission lines which 
have been observed to be strongest in Seyfert 
galaxies are the recombination lines of hydrogen, 
Pa/3 (A=1.2818/zm) and Br 7 (A=2.1655/zm), col- 
lisionally excited lines of singly ionized iron, e.g., 
[Fe II](A=1.2567//m), and rotational- vibrational 
lines of molecular hydrogen, e.g., H 2 ^=1-0S(1) 
(A=2.1212/im). Additionally, there are very high 
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excitation coronal lines, e.g., [S IX] (A=1.2525/xm), 
which have been seen in a number of Seyfert galax- 
ies (e.g., Moorwood & Oliva 1994; Knop et al. 
1996). To date, most published infrared spectra 
of Seyfert 2 galaxies use single large beams cen- 
tered on the nucleus (e.g., Goodrich et al. 1994). 
Thus, very little spatial information on the strong 
near-infrared lines is available. However, some 
studies have shown that near-infrared line emis- 
sion is extended over a spatial scale easily probed 
with ground-based spectroscopy. For example, 
near-infrared [Fe II] and H2 emission has been 
seen extended over several arcseconds (<1 kpc) in 
the nucleus of the Seyfert 2 galaxy NGC 1068 (Bli- 
etz et al. 1994), the Seyfert 1.5 galaxy NGC 7469 
(Genzel et al. 1995), and the Seyfert 1.5 galaxy 
NGC 4151 (Knop et al. 1996). In each case, there 
were overall spatial and kinematic similarities be- 
tween the infrared and optical line emission, but 
differences were also apparent on spatial and ve- 
locity scales within the NLR. Clearly, there is 
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important diagnostic information about the cir- 
cumnuclear regions of Seyfert 2 galaxies to be 
gained from spatially resolved near-infrared spec- 
troscopy 

This paper presents longslit near-infrared spec- 
tra of the Seyfert 2 galaxies Mk 1066, NGC 2110, 
NGC 4388, and Mk 3, in wavelength ranges which 
include the lines Pa/3, Br7, [Fe II] (A=1.2567/xm), 
H 2 z/=l-0 S(l), and [S IX](A=1.2525/xm). These 
four objects were chosen from a larger project to 
search for spatially extended infrared line emission 
in ^20 Seyfert 1 and Seyfert 2 galaxies. Because 
desirable targets in this larger project were those 
likely to show extended emission, Seyfert galax- 
ies visible from Palomar Observatory (5 > —20°) 
were selected with redshift cz<5400 km s _1 , corre- 
sponding to a spatial scale of 350 pc/". (Through- 
out this paper, for the determination of distances 
and hence projected scales, we assume Ho = 75 km 
s _1 Mpc -1 .) Galaxies with known extended opti- 
cal ionized gas were preferentially observed, fol- 
lowed in priority by galaxies with spatially ex- 
tended radio emission on scales of > 1". The final 
sample of galaxies observed includes 1 1 Seyfert 1 .x 
(1, 1.5, and 1.9) galaxies and 12 Seyfert 2 galax- 
ies. This paper reports the results from this pro- 
gram for the four Seyfert 2 galaxies that show the 
best well-detected extended circumnuclear emis- 
sion line regions in the near-infrared. 

Throughout this paper, for the determination 
of distances and hence projected scales, we assume 
H = 75 km s -1 Mpc" 1 . 

2. Observations and Data Reduction 

The near-infrared spectra discussed in this pa- 
per were obtained at the Hale 200" telescope with 
the Palomar Near-Infrared Spectrometer (Larkin 
et al. 1996). Observations at resolution A/ A A ~ 
1200 were made in two spectral ranges for each 
galaxy. One spectral range in the J-band atmo- 
spheric window covers 0.07/im and includes the 
[S IX], [Fe II], and Pa/3 lines in the rest frame of 
each galaxy. The other spectral range in the K- 
band atmospheric window covers 0.13/xm and in- 
cludes the H 2 i/=l-0 S(l) and Br7 lines. Table 1 is 
a log of the observations, including the resolution 
(in km s -1 ), on-source integration time, and posi- 
tion angle of the slit for each observation. Since we 
are interested in exploring the near infrared emis- 



sion line flux ratios and dynamics as a function 
of distance from the Seyfert nucleus, we chose slit 
position angles to coincide with the position an- 
gle of known ionization cones or other anisotropy 
previously mapped through narrow band optical 
emission line imaging. In the case of two of the 
galaxies (Mk 1066 and NGC 4388), there are addi- 
tional observations along a secondary slit through 
the nucleus, perpendicular to the primary position 
angle. 

The spectrometer is a reimaging system which 
disperses the light from the telescope and feeds the 
spectra to an externally mounted infrared cam- 
era. The spatial scale of the 256x256 NICMOS 
HgCdTe array in the camera is 0.165" pixel -1 . 
The width of the adjustable slit was approximately 
0.6", yielding the resolutions listed in Table 1. 
In order to accurately center the nucleus in the 
slit, a fold mirror inside the spectrometer, but in 
front of the grating, allows infrared imaging with a 
10" x 40" field of view. Between subsequent obser- 
vations of each galaxy, the nucleus was moved back 
and forth along the slit by 20", to allow for back- 
ground subtraction while maximizing the amount 
of on-source integration time. Each individual ex- 
posure was 600s for the J-band (1.2/im) wave- 
length range, and 300s for the K-band (2.2/mi) 
wavelength range. Along with the observation of 
each galaxy, we observed a G-dwarf star at a sim- 
ilar airmass (within 0.1 of that of the galaxy), and 
where possible in the same area of the sky, in order 
to correct for atmospheric transmission. A chop- 
ping secondary mirror on the telescope uniformly 
spread the light of this star along the length of the 
slit. In this way, the spectrum of the atmospheric 
calibrator is also used as a flat-field. 
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Galaxy 


Wavelength 


On-Source 


Slit 


Res. 


psf 


Date of 






Integration 


PA 


km s 1 


r VV rllvl 


Observation 


Mk 1066 


1.25/im-1.31/nm 


2400s 


135° 


242±15 


0.7 


1995/09/09 


Mk 1U66 




z400s 


/ICO 

40 


000_l_ 1 EC 


0.6 


i nnfi /an /in 
1996/09/ 10 


Mk 1066 


2.10^m-2.23/um 


2400s 


135° 


289±19 


0.6" 


1995/09/08 


Mk 1066 




2400s 


45° 


289±19 


0.5" 


1995/09/10 


NGC 2110 


1.25/im-1.31/im 


2400s 


160° 


229±16 


~ 1.1" 


1996/01/05 


NGC 2110 


2.09Mm-2.22/im 


1800s 


160° 


285±19 


~ 1.1" 


1996/01/05 


NGC 4388 


1.25/xm-1.31^m 


2400s 


30° 


229±15 


~ 0.8" 


1996/01/03 


NGC 4388 




2400s 


120° 


244±15 


0.8" 


1996/04/08 


NGC 4388 


2.10/wn-2.23//m 


2400s 


30° 


261±19 


~ 0.7" 


1996/01/03 


NGC 4388 




1800s 


120° 


303±19 


0.8" 


1996/04/08 


Mk 3 


1.26pim-1.32/im 


3000s 


113° 


243±15 


~ 1.0" 


1995/11/06 


Mk 3 


2.14^m-2.27/um 


1800s 


113° 


292±19 


~ 1.0" 


1995/11/07 



Table 1: Observation log. References for position angles are: Mk 1066: Bower et al. (1995); NGC 2110: 
Mulchaey et al. (1994); NGC 4388: Corbin, Baldwin, & Wilson (1988); Mk 3: Mulchaey et al. (1996). 
The resolution column indicates the FWHM of an unresolved line measured from OH sky lines in one of 
the galaxy frames. The point spread function (psf) of the atmospheric seeing was determined for Mk 1066 
and for the April 1996 observations of NGC 4388 by observing a star, guided in the same manner as the 
spectra, in between spectral observations. For Mk 3, the seeing was estimated based on 2.2/xm imaging of a 
quasar somewhat later than the observations of Mk 3. For NGC 2110 and the January 1996 observations of 
NGC 4388, unguided G-star and photometric standard star images were used to estimate the seeing. 
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The data were reduced following the procedure 
described in Knop (1997). Bad pixels were re- 
moved by linear interpolation, and the spatial and 
spectral curvature of the field of view was mea- 
sured and removed using a point source spectrum 
and the spectrally unresolved OH night sky lines. 
Spectral data frames were divided by the appro- 
priate G-star spectrum, in order to remove at- 
mospheric absorption lines and to correct for effi- 
ciency variations across the field of view. In order 
to correct for the blackbody slope of the G-star 
continuum, and to remove absorption features in- 
trinsic to the star, each G-star spectrum was first 
divided by a solar spectrum following the proce- 
dure of Maiolino, Rieke, & Ricke (1996). 

3. Results and Discussion: Individual Ob- 
jects 

3.1. Mk 1066 

Mk 1066 is an inclined SB0+ galaxy with a 
systemic velocity of 3625 km s _1 (Bower et al. 
1995). This corresponds to a calculated dis- 
tance of 48 Mpc, which yields a spatial scale of 
234 pc arcsec -1 . Hubble Space Telescope (HST) 
imaging of the nuclear regions shows a narrow 
"jetlike" feature in a narrow band image which 
includes [OIII] and H/3. This jet is oriented at 
315° and extends ~ 1.4" northwest of the nucleus. 
In Ha and [Nil], the jet is observed on both sides 
of the nucleus (Bower et al. 1995). The nuclear 
radio source is linear, and extended over 2.6". The 

5 GHz radio emission is oriented at the same po- 
sition angle as the optical line emission (Ulvestad 

6 Wilson 1989). 

We obtained infrared spectra of Mk 1066 along 
two position angles, 135/315° and 45/225°. The 
135° angle of the slit was chosen to align with 
the position angle of the optical ionization cones 
(Bower et al. 1995). These spectra show strong 
line emission, concentrated on the nucleus. Four 
emission lines are visible in Figures 1, 2, and 3: 
Pa/3 and [Fe II] in the J-band spectra, and Br7 and 
H 2 i/=l-fl S(l) in the K-band spectra. No [S IX] 
(A=1.2525/im) is visible on the nucleus (to an up- 
per limit on the equivalent width of 2 x 10~ 4 /im). 
In addition, there is data along a 45° slit chosen 
to be perpendicular to the primary slit. The in- 
frared lines are visible over a spatial extent of <~ 5" 
(1.2 kpc) along the 135° slit, in comparison to 



the FWHM of the nuclear continuum which was 
<~ 1.8" at 1.25/im and <~ 1.1" at 2.15^m in our 
spectra (both of which are greater than the 0.5- 
0.7" FWHM of the atmospheric seeing). In con- 
trast, the lines are less spatially extended along 
the 45° slit, and are visible only over an extent of 
< 3". The Br7 line is strongest on the nucleus, 
but the Pa/3, [Fe II], and H 2 lines peak about 0.5" 
northwest of the nucleus (see Figure 5). 

Figure 1 shows two greyscale plots of the re- 
duced 1.25-1.31^m spectra along both position 
angles. The horizontal axis is observed wave- 
length, and the vertical axis is position along the 
slit in arcseconds from the nucleus. The Pa/3 and 
[Fe II] lines are spatially resolved along both po- 
sition angles in Figure 1. The tilt in the lines in- 
dicates that there is a spatial velocity gradient in 
the [Fe II] and Pa/3 lines along the 135° slit; this 
will be discussed in section 3.1.2. Along the 45° 
slit, no velocity gradient is immediately apparent. 
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Fig. 1.— Spectra of Mk 1066 in the 1.25 - 1.3 l^m 
range. Position is distance from the nucleus in 
arcseconds along the indicated position angle of 
the slit. 

Figures 2 and 3 are one-dimensional projec- 
tions in various spatial bins extracted from the 
two-dimensional data. Figure 2 shows the spec- 
tra along the 135° slit, and Figure 3 along the 45° 
slit. The size of each spatial bin used in the ex- 
tractions is 4 pixels, or 0.67". This matches the 
point spread function (psf) at the time of the ob- 
servations (0.5"-0.7"). 

The Pa/3 profile near the nucleus shows en- 
hanced emission blueward of the line center. These 
blueward features on Pa/3 are visible not only in 
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Mk 1066, but also in several other Seyfert galax- 
ies (Knop et al. 2000). Were this Pa/3 emission, 
it would be at a velocity of ~ —700 km s _1 rela- 
tive to the centroid of the Pa/3 emission, which is 
higher than the velocities observed in other lines 
in Mk 1066. These blue features on Pa/3 are prob- 
ably primarily due to a blend of two He I lines 
at 1.2785/mi and 1.2791/im, which have been re- 
solved from Pa/3 in high resolution spectroscopy 
of the planetary nebula BD+30°3639 (Goodrich 
et al. 1994). There may also be a small contri- 
bution from [Fe II] (A = 1.2788//m). For T c ~ 
10 4 K and n e ~ 10 4 cm' 3 , the strength of the 
[Fe II] (A = 1.2788/xm) line is expected to be 
< 0.1 of the strength of the [Fe II] (A=1.2567/zm) 
line (Goodrich et al. 1994). In Mk 1066, [Fe II] 
(A = 1.2788/im) could account for at most half of 
the emission on the blue side of Pa/3. 

3.1.1. Integrated Line Flux Ratios 

Table 2 in Appendix A lists the equivalent 
widths and flux ratios of the lines observed in 
Mk 1066. Because [Fe II] and Pa/3 were observed 
simultaneously and are close in wavelength, we 
may obtain accurate [Fe II] /Pa/3 flux ratios which 
are relatively insensitive to uncertainties in flux 
calibration or differential reddening. Similarly, we 
may obtain reliable H 2 /Br7 flux ratios from our 
K-band spectra. Figure 4 shows these ratios as a 
function of position along the slit, using the total 
integrated fluxes for each species in each spatial 
bin to derive the line flux ratios. 

A striking feature of Figures 3 and 4 is the 
rapid rise of the H 2 /Br7 flux ratio with distance 
from the nucleus along the 45° slit, perpendicu- 
lar to the position angle of the optical ionization 
cones. The H 2 line is detected 2" away from the 
nucleus in both directions, while the Br7 line is 
detected only within 0.67" of the nucleus. The 
H 2 /Br7 flux ratio 1.33" from the nucleus is >3, in 
contrast to the nuclear value of 0.8. By compari- 
son, along the 135° slit, the H 2 /Br7 flux ratio al- 
ways remains between about 0.4 and 0.9, with an 
indication of a gradient across the nucleus. The 
[Fe II] /Pa/3 flux ratio varies with position along 
the 45° slit, although not as rapidly as does the 
H 2 /Br7 flux ratio. There is a clear local minmum 
in the [Fe II] /Pa/3 ratio at the nucleus along both 
slits. 
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Fig. 2.— Spectra of Mk 1066 along the 135° slit. 
Each spectrum represents the spectrum in a rect- 
angular beam which is 0.67" (4 pixels) along the 
slit by the width of the slit (~ 0.6"). Spatial bins 
are adjacent. The position of the center of each 
spatial bin, relative to the nucleus, is indicated in 
the lower left of each panel. Flux units (Fa) are 
arbitrary. 
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Fig. 3.— Spectra of Mk 1066 along the 45° slit. 
Each spectrum represents the spectrum in a rect- 
angular beam which is 0.67" along the slit by the 
width of the slit <~ 0.6". Spatial bins are adjacent. 
Flux units (F\) are arbitrary. 



3.1.2. Velocity structure 

Figure 5 is a set of position versus velocity plots 
along the 135° slit for each of the four lines identi- 
fied in Figures 2 and 3. In each smoothed 
continuum has been subtracted from the spec- 
trum. Velocity is plotted relative to 3625 km s -1 , 
the systemic velocity of Mk 1066 (Bower et al. 
1995), and position is distance in arcseconds on 
the sky from the continuum peak. A velocity gra- 
dient is visible in each of the lines, in the sense 
that the lines are blucshifted southeast of the nu- 
cleus and rcdshiftcd northwest of the nucleus. The 
slope of the velocity gradient does not appear 
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Fig. 4.— Flux ratios for Mk 1066. The two plots 
on the left are for the 135° slit, the two plots on 
the right for the 45° slit. 



to be the same for all of the lines. Specifically, 
the H 2 line appears to have a steeper gradient 
(~300 km s -1 arcsec" 1 , or ^1300 km s" 1 kpc" 1 ) 
in comparison to the other three lines (whose 



gradients are ~200 km s" 1 arcsec -1 , or ^900 
km s _1 kpc" 1 ). In addition, Br7 appears to have a 
slightly steeper gradient (by ^50 km s -1 arcsec -1 ) 
than Pa/3 southeast of the nucleus. The gradients 
of the [Fe II] and Pa/3 lines are similar. 

Another feature apparent in Figure 5 is that all 
of the lines, except for Br7, peak northwest of the 
nucleus. This is discussed in section 3.1.4. 

Table 6 in Appendix A lists the velocity cen- 
troids and full widths at half maxima of the com- 
plete line profiles for the lines observed in Mk 1066 
along both slits in each spatial bin. Many of 
the lines show a complicated asymmetric profile. 
These lines were decomposed into the smallest 
number of Gaussian components necessary to sta- 
tistically describe the line profile. Figure 6 shows 
plots of the velocity of the strongest narrow Gaus- 
sian component of each line as a function of posi- 
tion along the slit. Usually, this is equivalent to 
a plot of the mode (peak) of the line profile as a 
function of position, because most of the line pro- 
files in Mk 1066 are dominated by a strong, central 
narrow component. For reference, a line connect- 
ing the Br7 velocity peaks is drawn on each plot, 
and the Ha velocity curve (Bower et al. 1995) is 
shown on the Pa/3 plot. 
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Fig. 6. — Velocity of the primary narrow compo- 
nent of each line versus position along the 135° 
slit for Mk 1066. The data for Br7 are shown as 
a dotted line in each plot. The dashed line on the 
Pa/3 plot is the Ha data of Bower et al. (1995). 



Fig. 5. — Position versus velocity plots for 
Mk 1066. Position along a 135° slit. Velocity is 
relative to 3625 km s _1 , the systemic velocity of 
Mk 1066 (Bower et al. 1995). 



The same relative velocity gradients between 
the different species mentioned for Figure 5 are 
visible in Figure 6. It is apparent the difference 
between the H 2 and Br7 velocity slopes is greater 
than the difference between the H 2 and Pa/3 ve- 
locity slopes. Two arcseconds southeast of the nu- 
cleus, all of the velocity profiles except for that of 
Br7 (which is undetected) show a local velocity 
reversal of 50-75 km s^ 1 . 

In addition to the primary narrow emission line 
components whose centroid velocities are plotted 
in Figure 6, there are weaker components which 
are displaced in velocity relative to the center of 
the lines. Southeast of the nucleus, where the 
peaks of the line profiles are blucshiftcd relative 
to the systemic velocity, all of the lines show a 
noticeable enhanced red wing (although it is very 
weak in H 2 ). On the nucleus and to the northwest 
where the peaks of the line profiles are rcdshiftcd, 
the Pa/3 and [Fe II] lines show an enhanced blue 
wing. While this blue wing is almost certainly a 
feature of the [Fe II] line profile itself, the blue 
wing of the Pa/3 line suffers contamination from 
other species as discussed above. 



For illustration, Figure 7 shows the profiles of 
the [Fe II], Pa/3, and Br7 lines on the nucleus 
and 1.33" southeast and northwest of the nucleus. 
Southeast of the nucleus, the red wing is stronger 
relative to the rest of the line in [Fe II] and Br7 
than it is in Pa/3. Similarly, northwest of the nu- 
cleus, the blue wing is most prominent in [Fe II]. 
It is not visible in the Br7 line, and in Pa/3 it may 
be contaminated by the aforementioned He I and 
[Fe II] lines. 

3.1.3. Comparison with Optical Spectroscopy 

Bower et al. (1995) present longslit optical spec- 
troscopy of Mk 1066 along a 134° slit. The Ha, 
H/3, [Nil], [SII], and [OI] lines all have a velocity 
gradient extending from +100 km s _1 2-3" north- 
west of the nucleus to -90 km s _1 2-3" southeast of 
the nucleus. This is similar to the velocity struc- 
ture seen in the infrared lines (see Figure 6), al- 
though the peaks of the Pa/3 and particularly the 
Br7 lines are slightly blueshifted relative to the 
optical lines to the southeast of the nucleus. 

In addition to the primary kinematic compo- 
nent seen in optical emission lines, Bower et al. 
(1995) identify a second kinematical component, 
most clearly seen in [OIII]. The [OIII] emission 
is slightly blucshiftcd relative to the other lines 
northwest of the nucleus and redshifted by as much 
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Fig. 7. — Selected line profiles of Mk 1066 along 
the 135° slit, chosen to illustrate shoulders seen in 
asymmetrical line profiles. One resolution clement 
is approximately 3 bins on the plot. 



as 200 km s _1 relative to the other lines southeast 
of the nucleus. This structure is qualitatively sim- 
ilar to the enhanced wings seen in the infrared 
lines discussed in section 3.1.2. Specifically, the 
infrared lines show an enhancement blueward of 
the peak on the nucleus and to the northwest, and 
they show an enhancement redward of the peak to 
the southeast. It is likely that some of the wings 
seen on the infrared emission lines are associated 
with the kinematic component Bower et al. (1995) 
identify in [OIII]. 

3.1.4- Extinction 

Since Pa/3 and Br7 are both hydrogen recombi- 
nation lines (n=5-3 and n=7-4 respectively), they 
should arise in the same physical locations. Their 
intrinsic Case B flux ratio is well known (assum- 
ing radiative transfer effects are negligible). Be- 
cause of reddening, in dusty regions the observed 
Pa/3 line will be suppressed relative to the observed 
Br7 line. If the dust is intermixed with the line 
emitting clouds, and there are strong velocity gra- 
dients, the two lines can appear to have different 
velocity structures. In Mk 1066, both the infrared 
and optical data indicate that extinction effects 



are stronger to the southeast than they are to the 
northwest. 

A difference in both the spatial and velocity dis- 
tributions of Pa/3 and Br7 is visible in Figure 5. 
The Br7 line is peaked close to the continuum 
peak, and is stronger relative to Pa/3 toward the 
southeast than it is toward the northwest. The 
simplest explanation of this difference is a larger 
extinction in the southeast. Spatially variable ex- 
tinction may also be responsible for the slightly 
different Pa/3 and Br7 velocity gradients seen in 
Figure 5. To the southeast, where dust extinction 
appears greater, the discrepancy between the ve- 
locity gradients for the Pa/3 and Br7 lines is larger. 
Extinction could also explain the difference be- 
tween the velocities of the infrared lines and the 
optical Ha line southeast of the nucleus (Bower et 
al. 1995). The Ha- infrared slope discrepancy is in 
the same sense as the Pa/3-Br7 slope discrepancy, 
i.e., lines at shorter wavelengths have shallower 
gradients; see Figure 6. 

This same extinction effect is seen in the red 
wings of the Pa/3 and Br7 lines southeast of the 
nucleus as discussed in section 3.1.2; the wing is 
relatively more prominent in Br7 than in Pa/3 sug- 
gesting that the extinction to the region where the 
wings originate is greater than the extinction to 
the region where the bulk of the line emission orig- 
inates. If the wings represent outflowing gas, then 
the redshiftcd receding gas, which may be par- 
tially obscured by dust in the disk of the galaxy, 
would show higher extinction. There is also evi- 
dence from optical wavelengths that the outflow- 
ing gas in the southeast is heavily extinguished. 
The optical line emission is seen emerging from 
the nucleus in ionization cones at position angles 
of 135°/315°. The [OIII]+H/3 narrowband image 
(Bower et al. 1995) is much stronger to the north- 
west, suggesting suppressed emission to the south- 
east. The 5 GHz radio emission (Ulvestad & Wil- 
son 1989) indicates that the strength of the out- 
flow is slightly asymmetric, but the difference in 
the strength of the radio emission northwest and 
southeast of the nucleus is not as great as the dif- 
ference in the strength of the [OIII] emission which 
Bower et al. (1995) associate with the outflow. 

Extinction cannot explain the discrepancies be- 
tween the H2 and Br7 lines, since the differential 
extinction between the two is negligible. The H 2 
line shows a different velocity gradient than Br7 
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(Figure 5), and the peak of the H2 emission is 
shifted to the northwest relative to the peak of 
the Br7 emission. The differences in the structure 
of these lines may indicate that the emission seen 
from each species originates from different clouds 
along the same line of sight. Additionally, there 
may be a systematic variation in the physical con- 
ditions along the cone axis of the galaxy. 

3.1.5. A Picture of the Circumnuclear Regions of 
Mk 1066 

The infrared data in combination with the data 
in Bower et al. (1995) indicate that there are two 
dominant kinematic components in near the nu- 
cleus of Mk 1066. One component shows up in 
infrared H 2 , in the peaks of the infrared lines, and 
in optical Ha, [N II], and other species. This com- 
ponent shows a clear velocity gradient along the 
135° slit, is blueshiftcd to the southeast, and is 
redshifted to the northwest. The second, higher- 
excitation component shows up in the wings on 
the infrared lines and in optical [O III], and has 
the opposite velocity trend from the first compo- 
nent. 

Bower et al. (1995) identify the first, low- 
excitation component with a rotating disk. The 
velocity gradient seen in the infrared line emission 
for this component along the 135° slit is consis- 
tent with a rotation. The much smaller (though 
non-zero) velocity gradient along the 45° slit in- 
dicates that the major axis of the rotating disk 
must be closer to 135° than 45°; a fit of a uni- 
formly rotating disk to the peaks of the infrared 
line emission along both slits suggests that the 
major axis of the disk is at ~ 120°. Following 
Bower et al. (1995), we identify the second com- 
ponent with high-excitation outflowing gas. This 
second component is not seen off of the nucleus 
in the infrared spectra along the 45° slit, and the 
optical [O III] emission is clearly extended along 
an angle of ~ 135°. Consequently, the data indi- 
cate that the outflow is not perpendicular to the 
rotating disk; rather, the two physical systems are 
at low relative inclination when projected on to 
the plane the sky. 

Figure 8 is a cartoon drawing of the system 
based on the optical and infrared data. In this 
picture, there is an inclined disk with a major axis 
of 120° projected on the sky. This orientation 
explains both the greater spatial extent and the 




Fig. 8. — Cartoon drawing of the nuclear regions 
of Mk 1066, from the point of view of the ob- 
server. The rotating disk is the source of low ex- 
citation emission, and the source of the primary 
velocity gradient seen in most optical and infrared 
lines. The outflow is higher excitation emission, 
and shows up as wings in the infrared lines. The 
dusty disk extinguishes features from the south- 
east cone of the outflow. Superimposed on the 
cartoon are the angles of the 135° and 45° infrared 
slits. 



larger velocity gradient in the infrared lines along 
the 135° slit. The disk is dusty, explaining the ex- 
tinction effects discussed in section 3.1.4. There 
is a biconical outflow associated with the higher- 
excitation gas, as drawn. To the northwest, this 
outflow is relatively unobscured, as it is in front of 
the plane of the galaxy. The blue wings seen on 
the infrared lines to the northwest are associated 
with this outflow. To the southeast, the outflow is 
obscured by the plane of the galaxy, suppressing 
the [OIII] emission observed in that direction. Ex- 
tinction resulting from this obscuration enhances 
the Br7/Pa/3 ratio on the red side of the line pro- 
file southeast of the nucleus (see Figure 7). 

3.1.6. [Fe II] and H 2 Emission Processes 

The fact that various features of the [Fe II] 
emission line profiles correspond to those identi- 
fied with the active nucleus in other species, such 
as the wings seen in [OIII] (Bower et al. 1995), 
suggests that the [Fe II] emission is dominated by 
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processes directly associated with the Seyfert nu- 
cleus. The radio emission of Mk 1066 (Bower et 
al. 1995) appears to represent two oppositely di- 
rected radio jets aligned with the outflow seen in 
[OIII]. In the wings on the infrared line profiles 
we associate with this outflow, [Fe II] emission is 
enhanced relative to Pa/3 (see Figure 7). Shocks 
resulting from the interaction of outflowing radio 
plasma with the ambient gas are a natural mech- 
anism for the production of [Fe II] emission. Fast 
shocks can both destroy dust grains, increasing 
the gas phase abundance of iron, and create large 
partially ionized regions, two conditions which en- 
hance [Fe II] emission. 

The most striking feature of the H2 emission 
is that its strength does not drop off as quickly 
with distance as the other infrared lines from the 
nucleus perpendicular to the outflow axis. While 
Br7 is only detected out to 150pc from the nu- 
cleus at a position angle of 45°, H 2 is detected to 
>300pc. The presence of H 2 emisison along both 
slits shows that warm molecular gas is ubiquitous 
in the near-nuclear disk of Mk 1066. Since the H 2 
is not confined to the cone axis, it is evidently not 
directly dependent on ionization from the active 
nucleus or the radio jets for its excitation. There 
is probably little H 2 emission associated with the 
high excitation outflow, since the wings seen in the 
[Fe II] and Pa/3 lines associated with this outflow 
were not detected in the H 2 line profiles. 

Heating by X-rays from the Seyfert nucleus 
likely cannot explain the extent of the H 2 emission 
perpendicular to the ionization cone axis. Slow 
(v< 25 km s _1 ) shocks originating in the nu- 
cleus may heat molecular clouds in the disk of the 
galaxy immediately surounding the nucleus, but 
the relatively flatter distribution of H 2 relative to 
the other species along the 45° slit argues for an 
extended source of heating for the molecular gas 
in the regions within ~ 500 pc of the nucleus of 
Mk 1066. 

3.1.7. Summary 

1. The infrared Pa/3, Br7, and [Fe II] lines are 
extended over 5" (1.2 kpc projected on the 
plane of the sky) along a 135° position angle 
coincident with the direction of the optically 
identified ionization cones. Along a perpen- 
dicular angle, their intensity drops off faster 
with distance from the nucleus; 1.3" (0.3 kpc 



projected) away from the nucleus at 45°, the 
intensity of the lines relative to the nucleus 
is 20% of their value along a position angle 
of 135°. 

2. The infrared lines have strongly peaked nar- 
row components. The peaks of the lines have 
a velocity gradient similar to the optical Ha, 
H/3, [Nil], [SII], and [OI] lines along a posi- 
tion angle of 135°, which Bower et al. (1995) 
identify as part of a rotating, LINER-like 
disk. 

3. There are secondary components in the form 
of enhanced wings in the infrared lines whose 
kinematics differ from the strong, narrow 
components in a manner similar to the differ- 
ence between the kinematics of [OIII] emis- 
sion and the other optical lines. The [OIII] 
emission line and the near infrared line wings 
are likely associated with a nuclear outflow 
which produces blueshifted emission to the 
northwest of the nucleus. 

4. The relative spatial distributions of Br7 and 
Pa/3 indicate that there is greater extinction 
to the southeast of the nucleus than there 
is to the northwest. This may explain differ- 
ences seen in the velocity centers of hydrogen 
recombination lines to the southeast, which 
is most likely due to obscuration of the gas 
by an inclined disk in Mk 1066. 

5. The H 2 line is more extended than the other 
lines along an angle perpendicular to the ion- 
ization cone. Along the axis of the cone, it 
shows a steeper velocity gradient than Pa/3, 
Br7 or [Fe II]. 

6. We propose that there is a measurable con- 
tribution to the [Fe II] emission from shocks 
associated with the nuclear outflow. It may 
have a direct link to the radio jets seen at 
3.6cm. The H 2 emission, on the other hand, 
arises primarily from the inner (500pc) disk 
of Mk 1066, most of which is presumably 
shielded from direct illumination of the nu- 
cleus. This warm molecular gas may be 
heated by extended local sources such as star 
formation. 
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3.2. NGC 2110 



NGC 2110 is an SO galaxy with a systemic ve- 
locity of 2342 km s -1 (de Vaucouleurs et al. 1991), 
corresponding to a calculated distance of 31 Mpc, 
implying a spatial scale of 150 pc arcsec" 1 . HST 
narrowband imaging in optical Ha+[NII] and 
[OIII] emission shows a narrow jctlikc feature near 
the nucleus extending to the north at a position 
angle of ^340°, and a weaker opposing feature 
extending to the south at ^160° (Mulchaey et 
al. 1994). The feature curves into an "S" shape 
4" away from the nucleus in each direction. At 
1.49 GHz, NGC 2110 shows a triple radio source 
extended over 4" with a position angle of ^0° 
(Ulvestad & Wilson 1984), slightly different from 
the position angle of the optical line emission. The 
curvature of the "S" of the optical line emission 
is towards the direction of the radio emission. In 
X-rays, NGC 2110 has a 2-10 kcV luminosity of 
1.4 xlO 43 erg s _1 with a comparatively high ab- 
sorbing column density of 7.4 x 10 22 cm 22 (Rcichcrt 
et al. 1985). This X-ray luminosity is closer to that 
typical for a moderate luminosity Seyfert 1 galaxy 



than for a Seyfert 2 galaxy (~ 5x 10 erg s 



- 1 ) 



(Kriss et al. 1980). ROSAT High Resolution Im- 
ager (HRI) data show a spatially extended X-ray 
source, with secondary X-ray emission 4" north of 
the nucleus which may be associated with the soft 
X-ray excess seen by Weaver et al. (1995). 

Our longslit spectra of NGC 2110 along the 
ionization cones at a position angle of 160° (Fig- 
ure 9) show strong emission lines in [Fe II], Pa/3, 
H 2 u=l-0 S(l), and Br 7 . No [S IX] (A=1.2525^m) 
is detected on the nucleus (with an upper limit on 
the equivalent width of 1 x 10~Vm). Figure 9 
shows the summed spectra in spatially adjacent 
1" bins along the slit. These lines are visible on 
the nucleus, and are seen to be extended by > 6" 
(900 pc) . In comparison, the FWHM of the contin- 
uum is 2.4" at 1.2/xm and 1.5" at 2.2/xm, while the 
seeing at the time of the observations was ~ 1.1". 
The most striking feature of the NGC 2110 J-band 
spectrum is the strength of the [Fe II] line. On 
the nucleus, the [Fe II] /Pa/3 flux ratio is 8.1±0.8, 
which is greater than that observed in any other 
Seyfert galaxy in this survey (Knop et al. 2000). 
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Fig. 9.— Spectra of NGC 2110 along the 160° 
slit. Each spectrum is produced from a 1.0 x 0.6" 
rectangular beam. Spatial bins are adjacent. 
The distance of each spatial bin from the nu- 
cleus in a southerly direction is indicated in each 
panel. Note the strength of [Fe II] and the large 
[Fe II] /Pa/3 flux ratio on the nucleus. Flux units 
(Fa) are arbitrary. 
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3.2.1. Integrated Line Flux Ratios 

Table 3 in Appendix A lists the equivalent 
widths and flux ratios of the lines observed in 
NGC 2110. The [Fe II] /Pa/3 and H 2 /Br 7 line ra- 
tios are very large on and near the nucleus. These 
flux ratios are plotted as a function of position 
along the slit in Figure 10. Near the nucleus, the 
H 2 /Br7 flux ratio is close to 3, locally rising at a 
position 2" north of the nucleus, and dropping off 
in the furthest spatial bins. The [Fe II] /Pa/3 flux 
ratio peaks on the nucleus, falling rapidly with ra- 
dius away from the nucleus. 
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Fig. 10. — NGC 2110 flux ratios as a function of 
position along the 160° slit. Positive positions are 
to the south. 



lines seen in NGC 2110. A smoothed continuum 
has been subtracted from the data in each case. 
Velocity is plotted relative to the systemic velocity 
of NGC 2110. 
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For most Seyfert galaxies, values of [Fe II] /Pa/3 
typically range between 0.5 and 2.0 (e.g., Goodrich 
et al. 1994; Simpson et al. 1996b), while the 
H2 i^=l—0 S(l)/Br7 flux ratio typically lies be- 
tween 0.5 and 1.5 (e.g., Moorwood & Oliva 1988). 
The flux ratios observed in NGC 2110 are signif- 
icantly higher than these, with the [Fe II] line 
in particular being unusually strong in compar- 
ison to the Pa/3 line. Indeed, the value of 8.1 
for the [Fe II] /Pa/3 flux ratio on the nucleus 
is higher than that observed in any published 
Seyfert galaxy spectrum (e.g., Goodrich et al. 
1994; Forbes and Ward 1993; Mouri et al. 1993). 
This high [Fe II] /Pa/3 ratio is due to both a strong 
[Fe II] line and a weak Pa/3 line in comparison 
to other Seyfert galaxies. The equivalent width 
of the [Fe II] line on the nucleus, 25xl0~ 3 ^m, is 
larger by a factor of 2-4 than that typically seen 
in Seyfert galaxies, although it is not the largest 
known equivalent width of [Fe II] in a Seyfert 2 
galaxy (e.g., NGC 449; Knop et al. 2000). 

3.2.2. Velocity structure 

In Figure 11, we show a set of position versus 
velocity plots for each of the four near-infrared 



Fig. 11. — Velocity versus position along 160° slit 
for NGC 2110. Velocity is relative to 2340 km s" 1 , 
the systemic velocity of NGC 2110. 

All of the lines in Figure 11 show a gradient 
with velocity increasing from north to south. Un- 
like Mk 1066, the velocity gradients appear to be 
quite similar for the four infrared lines seen in 
NGC 2110. Figure 12 is a plot of the velocity 
centroid of each line profile as a function of posi- 
tion for each of the four lines. The centroids for 
Pa/3 have excluded the flux due to the contami- 
nation from He I and [Fe II] on the blue side of 
Pa/3 (as discussed in section 3.1). Table 7 in Ap- 
pendix A lists these velocity centroids, as well a 
the full widths at half maximum of each of the 
line profiles. The velocity structure in Figure 12 
is dominated by the gradient from north to south. 
[Fe II] appears different from the other species in 
Figure 11 because it is much wider on the nucleus 
than it is off of the nucleus. The shape of the 
contours convolves the spatial variation of the line 
width and the gradient in the line's velocity. The 
slope of the gradient in all species is steeper to 
the south than to the north (Figure 12), reaching 
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velocities of about —100 to —150 km s 1 at 3". 
South of the nucleus and +250 to +300 km s _1 at 
3" South of the Nucleus. 
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Fig. 12. — Velocity centroid of the whole line 
profile versus position for each line at each po- 
sition along the 160° slit. Velocity is relative to 
2340 km s~\ the systemic velocity of NGC 2110. 



A second feature apparent in Figure 12 is that 
the [Fe II] line is systematically at a lower velocity 
than the other lines. This effect is strongest south 
of the nucleus, although here the strong blue wings 
in [Fe II] pull the centroids to lower velocities. In 
each spatial bin the velocity peak of the [Fe II] line 
is shifted to the blue with respect to the other lines 
by about 50-100 km s -1 , which is several times the 
uncertainty in the wavelength calibration. 

In many cases, the lines are asymmetric, al- 
though the asymmetries are mostly small enough 
that the velocity centroids are very close to the 
velocities of the peaks of the lines. On the nucleus 
and 1" south of the nucleus, Pa/3, [Fe II], and H 2 
show a blue wing in places, while north of the nu- 
cleus, Br7, [Fe II], and H 2 show red wings. 

The [Fe II] line is significantly broader than 
the other lines, particularly close to the nucleus 
(see Figure 11). On the nucleus, the FWHM of 
the [Fe II] line is 500 km s _1 , as compared to 
340 km s _1 for the Pa/3 line. In addition, the 
[Fe II] line has wings with FWHM ~ 850 km s" 1 
which are not visible in the Pa/3 line. The [Fe II] 
line is narrower at larger radii, having a FWHM 
of < 400 km s" 1 , 3" away from the nucleus. At 



3" south of the nucleus, the [Fe II] line remains 
< 100 km s _1 broader than the Pa/3 line, while 
north of the nucleus the line widths are compara- 
ble. 

3.2.3. Comparison with Previous IR Spectroscopy 

Storchi-Bergmann et al. (1999) (hereafter cited 
as S99) performed spatially extended infrared 
spectroscopy of NGC 2110 with resolutions of 
<~ 700 and <~ 2000. They have data along two 
slit angles: one at 170°, similar to our slit angle 
of 160°, and one at a perpendicular angle of 80°. 
The features of their spectra along the 170° slit 
are also seen in this work, including the notable 
increase in the [Fe II] /Pa/3 ratio on the nucleus. 
Any quantitative differences between the results 
of the two studies are easily explained by differ- 
ences in slit widths and seeings. One exception is 
the slope of the continuum in the 2.15^m spectra, 
where the slope seen in this work is opposite that 
seen in S99. This is probably also a result of dif- 
ferences in the sampled spatial region, but may be 
partially a result of differences in data reduction, 
indicating the difficulty in absolute calibration of 
the continuum in infrared spectra. Since this work 
is primarily concerned with the spatial and veloc- 
ity distribution of emission line radiation, it is 
important to consider how these differences affect 
calculated line flux ratios. For the differences seen 
between the two nuclear continuum spectra, the 
H 2 /Br7 flux ratio should not be misestimated by 
more than ^15%, comparable to our quoted er- 
ror on the nuclear flux ratio. S99 saw a nuclear 
H 2 /Br7 flux ratio of 3.0+2.0, consistent with our 
value of 3.2+0.6; this indicates that the flux ratios 
observed in each case are in fact robust. 

Along the perpendicular 80° slit, S99 have 
data in a wavelength range including Br7 and H 2 
(A=2.1212^m). H 2 is well detected to a greater 
radius along this 80° slit, which is perpendicular 
to the ionization cones. Although S99 have weak 
or no detection of Br7 off of the nucleus in cither 
direction, the S/N of the off- nuclear H 2 detection 
is better along the 80° slit than it is along the 
170° slit. This, coupled with our observations 
of the H 2 /Br7 ratio along our 160° slit, suggests 
that also as with Mk 3, the H 2 /Br7 ratio off of 
the nucleus is larger in the slit perpendicular to 
the ionization cone than it is along the cone axis. 
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3.2.4- Comparison With Optical Spectra and Ra- 
dio Morphology 

S99 report the detection of a broad Pa/? line in 
NGC 2110. This line has a width of 1200 km s _1 , 
and is highly asymmetric to the blue. Although 
the same spectral feature is present in the data 
here (see the nuclear spectrum in Figure 9), we 
interpret the enhanced blue wing as contamina- 
tion from He I and [Fe II] (see section 3.1), rather 
than as evidence for high velocity blueshifted Pa/3 
emission. The [Fe II] (A = 1.2788/xm) line which 
contributes to this contamination is generally only 
< 0.1 the strength of the [Fe II] (A=1.2567/zm) line 
(Goodrich et al. 1994). However, given that [Fe II] 
(A=1.2567/xm) is ~ 8 times as strong as Pa/3 on 
the nucleus of NGC 2110, it is likely that the [Fe II] 
(A = 1.2788/im) line contributes more significantly 
to the emission blueward of Pa/3 than is the case 
in most Seyfert 2 galaxies. 

The velocity structure seen in Ha, H/3, and 
[OIII] (Wilson et al. 1985) is quite similar to the 
velocity structure seen in the infrared lines plot- 
ted in Figures 11 and 12. Wilson et al. (1985) and 
Wilson & Baldwin (1985) describe this velocity 
field as a rotating disk with a projected major axis 
of 161°. (This paper only presents infrared data 
for NGC 2110 along one slit position angle. Al- 
though the rotation identified by Wilson & Bald- 
win (1985) is evident in this data, the lack of data 
along two spatial dimensions precludes any inde- 
pendent conclusions about the major axis of the 
rotation.) 

Near the nucleus, the [OIII] line has a higher 
root mean square (RMS) width (equivalent to la 
for a Gaussian line profile) than the H/3 line by 
100-200 km s" 1 (Wilson et al. 1985), although the 
higher [OIII] signal-to-noise allows weak wings to 
be visible in the former which are not visible in 
the latter. South of the nucleus, [OIII] shows a 
tendency toward enhanced blue wings, similar to 
what is observed in [Fe II] , and to a lesser extent in 
H2 . These enhanced wings show some spatial sim- 
ilarities with the 2 cm emission (Ulvestad & Wil- 
son 1984), which shows three peaks aligned along 
a roughly north-south axis. The strongest peak 
is on the nucleus. The blue wings are strongest 
in [Fe II], H 2 , and [OIII] 2" south of the nucleus, 
which is the approximate position of one of the 
secondary radio peaks. The other off-nuclear ra- 



dio peak is 2" north of the nucleus, where we see a 
red shoulder in [Fe II] and H 2 . This suggests that, 
as with Mk 1066, there may be two kinematical 
components: a rotating disk, and a second veloc- 
ity field, possibly an outflow associated with the 
radio jets (see below). Note that even though the 
projected position angle of these outflowing jets 
appears similar to the major axis of the rotation 
identified by Wilson & Baldwin (1985), the two 
physical systems need not be aligned in space. 

3.2.5. [Fe II] and H2 Emission Processes 

The fact that the [Fe II] /Pa/3 ratio is very large 
on the nucleus and rapidly drops with radius sug- 
gests that the [Fe II] emission is associated with 
the Seyfert nucleus itself. There is indirect evi- 
dence that the physical process responsible for the 
[Fe II] emission in NGC 2110 is X-ray heating and 
ionization. NGC 2110 has a 2-10 keV X-ray lu- 
minosity which is an order of magnitude greater 
than that typical for a Seyfert 2 galaxy. On the 
nucleus, the [Fe II] /Pa/3 ratio is also nearly an or- 
der of magnitude greater than the typical value for 
Seyfert 2 galaxies. Out to radii 3-4" north of the 
nucleus, the [Fe II] /Pa/3 ratio remains atypically 
large. About ~ 4" north of the nucleus, at the 
position of the secondary 0.1-2.4 keV X-ray emis- 
sion peak identified by Weaver et al. (1995), there 
is a lower limit on the [Fe II] /Pa/3 flux ratio of > 2 
(due to an undetected Pa/3 line). South of the nu- 
cleus, opposite the direction of the extended X-ray 
emission, the [Fe II] /Pa/3 ratio drops to unity. 

Although the [Fe II] excitation may be domi- 
nated by X-ray ionization, there is likely to be an 
additional contribution due to fast shocks result- 
ing from the interaction of outflowing radio plasma 
with surrounding gas. This would explain the en- 
hanced red and blue wings respectively north and 
south of the nucleus which are spatially coincident 
with the radio peaks seen by Ulvestad & Wilson 
(1984) at 1.49 GHz. 

Nuclear X-ray heating of molecular clouds may 
also be the primary mechanism responsible for the 
H 2 emission, although the enhanced H 2 emission 
is not as dramatic as is the case with [Fe II] in 
this galaxy. The H 2 /Br7 ratio in NGC 2110 is 
large for Seyfert galaxies, though the ratio is not 
strongly peaked on the nucleus as is the case for 
the [Fe II] /Pa/3 ratio. The H 2 /Br7 flux ratio is 
nearly constant over the 6" (900pc) where the lines 
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are detected. In contrast to Mk 1066, the presence 
of the wings in the H2 line profiles indicates that 
there may be warm molecular gas in the outflowing 
clouds. 

3.2.6. Summary 

1. In NGC 2110, the [Fe II] /Pa/3 flux ratio is 
unusually large, reaching a maximum value 
of 8.1 on the nucleus and dropping off rapidly 
with radius. The H 2 /Br7 ratio is also large, 
being > 3 within 2" of the nucleus. 

2. In general, the infrared lines show a smooth 
velocity gradient similar to the velocity gra- 
dient seen in optical emission by Wilson et 
al. (1985). 

3. The centroid of the infrared [Fe II] line is sys- 
tematically shifted with respect to the other 
lines by about 50-100 km s _1 to the blue. 
Additionally, the [Fe II] line is <~ 160 km s _1 
broader than the other infrared lines near 
the nucleus. 

4. There are enhanced wings on the infrared 
lines. To the south, where line centroids are 
redshifted with respect to the systemic ve- 
locity, there are enhanced blue wings. To the 
north, where line centroids are blueshiftcd 
with respect to the systemic velocity there 
are enhanced red wings. These wings are 
strongest in [Fe II] and H 2 at positions coin- 
cient with the secondary radio peaks of Ul- 
vestad & Wilson (1984). 

We suggest that the [Fe II] excitation is domi- 
nated by nuclear X-ray photoionization. However, 
the enhanced off-nuclear wings which are spatially 
coincident with the nonthermal 2cm radio peaks 
indicate that there is some emission associated 
with outflowing gas. Fast shocks may contribute 
to the [Fe II] emission at these locations. 

While the H 2 emission is likely to be enhanced 
by the strong X-ray flux in NGC 2110, the distri- 
bution of the warm molecular gas is not as cen- 
trally peaked as is the [Fe II] emission. This may 
indicate a large molecular reservoir surrounding 
the nucleus. In addition, the wings seen in the H2 
line profile suggest that at least some of the H 2 
emitting gas takes part in the outflow driven by 
radio jets. 



3.3. NGC 4388 

NGC 4388 is a highly inclined spiral galaxy at a 
systemic velocity of 2525 km s -1 (Petitjean & Dur- 
ret 1993). Although its Hubble type is uncertain 
(de Vaucouleurs et al. 1991), Corbin, Baldwin, & 
Wilson (1988) classify it as SB(s)b pec. It is gen- 
erally considered to be a member of the Virgo clus- 
ter (Hclou et al. 1980), implying that the distance 
to the galaxy is probably ~14 Mpc (Pierce & Tully 
1988). This yields a spatial scale of 68 pc arcsec -1 . 
Optical narrowband imaging (Corbin, Baldwin, & 
Wilson 1988) shows an Ha+[NII] extent of > 40" 
east and west of the nucleus, aligned with the ma- 
jor axis of the galaxy at position angle ~ 90°. 
This emission originates in the spiral arms of the 
galaxy, and is probably due to HII region com- 
plexes. Close to the nucleus on a spatial scale of 
~ 5", the [OIII] emission is extended about the nu- 
cleus in a halo-like distribution that has a major 
axis at a position angle of ~ 30°. High resolu- 
tion radio maps at 4.86 GHz show emission which 
is elongated southwest of the nucleus at a similar 
position angle of 200° (Hummel & Saikia 1991). 
At 4.86 GHz, NGC 4388 is double peaked, with 
a primary peak on the nucleus and a secondary 
peak 2.5" southwest of the nucleus. There is a 
weak tertiary peak about 1" northeast of the nu- 
cleus. NGC 4388 has also been detected in X-rays, 
with a 2-10 keV luminosity of <~ 2 x 10 42 erg s _1 
(Hanson et al. 1990). The soft (0.1-2.4 keV) X- 
ray emission of NGC 4388 is extended over > 15" 
(Matt et al. 1994). 

Because Ha emission is extended over > 40" 
along a position angle of 120° (Corbin, Baldwin, & 
Wilson 1988), we obtained a "sky" frame distant 
(200" east) from the center of the galaxy. Sub- 
tracting this sky frame from spectra of the galaxy 
verified that no infrared emission lines were de- 
tected at a radius greater than 5" in this direction 
to a 3<7 flux limit of <~ 2 x 10~ 16 erg cm~ 2 s _1 in 
a 1.0" x 0.6" spatial beam. 

Figures 13 and 14 are one-dimensional projec- 
tions in various spatial bins of the data along the 
slits oriented at position angles of 30° and 120° 
respectively. The size of each spatial bin used in 
the extractions is 6 pixels (1"), slightly larger than 
the psf at the time of the observations, which was 
0.7-0.8" for the spectra at both position angles. 
The nuclear spectra for the two slit position an- 
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gles show qualitative differences, which are likely 
attributed to vagaries of spatial sampling with the 
rectangular synthetic aperture. The 30° slit was 
chosen to align with the major axis of the distribu- 
tion of [OIII] emission near the nucleus (Corbin, 
Baldwin, & Wilson 1988). The 120° slit was cho- 
sen to be perpendicular to the [OIII] axis. 

In NGC 4388, the [Fe II], Pa/3, Br 7 , and H 2 
lines are detected with a spatial extent up to 8" 
(540 pc) along the 30° slit. Along the perpen- 
dicular 120° slit, the lines diminish quickly with 
distance from the nucleus, but are still visible out 
to a radius of 2-3". In addition to these lines, 
there is a narrow line seen on the nucleus blue- 
ward of [Fe II] , which we identify as the fine struc- 
ture line [S IX] at a rest wavelength of 1.25249/im 
(Kuhn, Penn, & Mann 1996). Oliva et al. (1994) 
identify this line as [S IX] in the Circinus galaxy. 
However, Simpson et al. (1996b) and Thompson 
(1995) identify this line as He I (A = 1.2528/xm), 
with possibly some contamination from [S IX] in 
NGC 1068 and NGC 4151 respectively. 

The arguments in favor of this line being pri- 
marily [S IX] are threefold. First, the He I line 
is not strong enough to account for the observed 
flux in the objects where this line is seen. Re- 
combination cascade calculations (Robbins 1968) 
indicate that the strength of He I (A = 1.2528/im) 
should be < 20% of the strength of the He I 
(A = 1.2785/mi) line. The latter line, blended with 
He I (A = 1.2791/mi) and [Fe II] (A = 1.2788/im), 
comprises emission which is frequently seen as a 
weak blue wing on the Pa/3 line profile. However, 
when seen (see also Knop et al. (2000)), the line 
blucward of [Fe II] is of comparable strength or 
stronger than the blend blueward of Pa/3, arguing 
against its identification as He I. 

Secondly, the observed wavelength centroid is 
more consistent with [S IX] than with He I. If 
the line is identified as [S IX], it is blueshifted 
by <~ 100 km s -1 relative to Pa/3. If the line is 
He I, its blueshift would be 174 km s _1 . Since 
the ionization potential of He is only 24.5 eV, we 
would expect He to show a similar spatial and ve- 
locity structure to the H+ lines. S 7+ , on the other 
hand, has a much larger ionization potential of 
328 eV (West et al. 1985), and may arise from 
very different physical regions than does most of 
the Pa/3 and Br 7 emission. 

Finally, the line we identify as [S IX] is only 



seen in a small number of galaxies for which in- 
frared spectra are available (Knop et al. 2000). A 
He I line would be expected to be seen with a 
similar strength relative to Pa/3 in most Seyfert 
galaxies, as the two lines should arise from similar 
physical regions. The strength of a high ioniza- 
tion coronal line may vary greatly between objects, 
however, due to variations in the ionization con- 
ditions, as well as differences in extinction to the 
regions responsible for the H + emission and the 
coronal emission. We therefore conclude that the 
line observed at a rest wavelength of 1.252^m is 
[S IX]. 

The [S IX] line is very strong on the nucleus. 
Its equivalent width of 4.5xl0~V m is nearly half 
that of [Fe II]. It is not as spatially extended as 
the other infrared lines along the position angle of 
either slit. It is marginally resolved spatially along 
the 30° slit, showing a stronger flux (0.39±0.08 of 
the nuclear [S IX] flux) in the ±1" spatial bins 
than would be expected for atmospheric seeing of 
0.8" (0.08 of the nuclear flux). It is not resolved 
along the 120° slit. 

3.3.1. Integrated Line Flux Ratios 

Table 4 in Appendix A lists the equivalent 
widths and flux ratios of the infrared lines detected 
in NGC 4388 along both slits. Figure 15 shows 
the [Fe II] /Pa/3 and H2/Br7 flux ratios along both 
slits. On the nucleus, the [Fe II] /Pa/3 flux ratio 
is about 0.4, a relatively low value for a Seyfert 
2 galaxy; however, the [Fe II] /Pa/3 flux ratio rises 
off the nucleus along both slits. The [Fe II] /Pa/3 
ratio is largest where the [Fe II] velocity centers 
are most discrepant from the Pa/3 velocity cen- 
ters, specifically to the southwest along the 30° 
slit and to the southeast along the 120° slit (see 
section 3.3.2). 

Along the 30° slit, which is coincident with the 
major axis of the optical [OIII] emission (Corbin, 
Baldwin, & Wilson 1988), the H 2 /Br7 ratio is 
<1 on the nucleus, decreasing slowly with radius. 
Along the 120° slit, however, the flux of the H 2 line 
does not diminish nearly as quickly as does the flux 
from the Br 7 line. At a distance 2" southeast and 
northwest of the nucleus, the Br 7 line is < 1/25 of 
its nuclear strength, while the H 2 line is 1/8 of its 
nuclear strength (see Figure 14). This behavior is 
qualitatively similar to that seen in Mk 1066. 
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Fig. 13.— Spectra of NGC 4388 along the 30° slit. 
Each spectrum is within a rectangular beam which 
is 1" along the slit by the width of the slit (<~ 0.6"). 
Spatial bins are adjacent. The distance of each 
spatial bin from the nucleus along a position angle 
of 30° is indicated in each spectrum. Flux units 
(Fa) are arbitrary. 
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Fig. 14.— Spectra of NGC 4388 along the 120° 
slit. Each spectrum is within a rectangular beam 
which is 1" along the slit by the width of the slit (~ 
0.6"). Spatial bins are adjacent. The distance of 
each spatial bin from the nucleus along a position 
angle of 120° is indicated in each spectrum. Flux 
units (F x ) are arbitrary. 



3.3.2. Velocity structure 

Figure 16 is set of continuum-subtracted posi- 
tion/velocity contour plots for the emission lines 
observed along the 30° slit for NGC 4388. There 
is an enhanced blue wing of Pa/3 on the nucleus, 
that is probably dominated by a combination of 
He I (A=1.2786 / um) and [Fe II] (A=1.2788/im), as 
discussed in section 3.1. 

Table 8 in Appendix A lists the centroid and 
width of each line in each spatial bin along both 
slits. Figures 17 and 18 show the velocity cen- 
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Fig. 15.— Flux ratios for NGC 4388. The two 
plots on the left are for the 30° slit; the two plots 
on the right, for the 120° slit. 



troids of the line profiles as a function of position 
along the 30° and 120° slits respectively. North- 
east of the nucleus along the 30° slit, the Pa/3, 
Br7, and [Fe II] lines agree in velocity centroid, 
and show only a small increase in velocity with 
distance from the nucleus. The H2 line is slightly 
redshifted relative to the other lines. Southwest 
of the nucleus, however, the lines show very differ- 
ent velocity structures. The Pa/3, Br7, and [Fe II] 
lines turn towards increasing velocity with radius 
to the southwest, while the velocity gradient of 
H2 flattens. The [Fe II] line curves back more 
sharply than the Pa/3 line, until its peak is about 
100 km s _1 redward of Pa/3 3" south of the nu- 
cleus. Figure 19 shows profiles of the Pa/3 and 
[Fe II] lines southwest of the nucleus, and indicates 
that the difference between [Fe II] and Pa/3 is not 
the result of errors in velocity measurements, but 
that the centroid and profile of [Fe II] is genuinely 
redshifted relative to Pa/3. Br7 also differs from 
Pa/3 southwest of the nucleus, increasing in veloc- 
ity closer to the nucleus than does Pa/3. 

Along the 120° slit, the infrared line fluxes drop 
off faster with radius than they do along the 30° 
slit. Along this slit, there is no evidence for signif- 
icant velocity changes in the Pa/3 line. Where it 
was detected, the velocity peak of the Br7 line is 
marginally consistent with that of the Pa/3 line, as 
is the [Fe II] line everywhere except for 1" south- 
east of the nucleus, where [Fe II] is blueshifted 
by 80 km s _1 relative to Pa/3. The H 2 line is al- 
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Fig. 16. — Position versus velocity plots for NGC 
4388. Position is along a 30° slit. Velocity is rela- 
tive to 2525 km s _1 , the systemic velocity of NGC 
4388. Contours in this figure are logarithmic and 
separated by a factor of two. 



ways blueshifted relative to the Br7 line by 50- 
100 km s _1 . 

In addition to different variations in velocity 
centroid with position, the lines have differing pro- 
files. Along the 30° slit, not only is [Fe II] gener- 
ally wider than Pa/3, and at a different velocity 
centroid southwest of the nucleus, but the asym- 
metries in the lines arc different. For instance, 
2-3" southwest of the nucleus the [Fe II] /Pa/3 ra- 
tio is preferentially enhanced at redshifted veloci- 
ties (Figure 19). These asymmetries suggest that 
systems of clouds with different flux ratios at dif- 
ferent velocities are superimposed along the line 
of sight. Along the 120° slit, the [Fe II] line has 
an enhanced red shoulder to the northwest of the 
nucleus. Again, the Pa/3 line is redshifted relative 
to the [Fe II] line, so this red shoulder may cor- 
respond to a component which is associated with 
the Pa/3 peak. 
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Fig. 19. — Line profiles of [Fe II] and Pa/3 in 
NGC 4388 2-3" southwest of the nucleus on the 
30° slit, and the log of the [Fc II] /Pa/3 ratio as a 
function of velocity. Points with greater than 50% 
uncertainty in the flux ratio have been omitted 
from the latter plot. 



Fig. 18. — Velocity centroid of each line versus 
position along the 120° slit for NGC 4388. Positive 
distance is to the southeast. 



3.3.3. Comparison with other observations 

In their nuclear spectrum which was marginally 
resolved spatially along a 13° slit, Winge et al. 
(2000) observed Pa/3 and [Fe II] with similar pro- 
files. This is in contrast to our data, where [Fe II] 
is clearly wider than Pa/? near the nucleus along 
a 30° slit. Narrowband optical imaging and spec- 
troscopy (e.g. Corbin, Baldwin, & Wilson 1988; 
Petitjean & Durret 1993), radio imaging (Hum- 
mel & Saikia 1991), as well as this study and 



that of Winge et al. indicate that NGC 4388 is a 
complicated system with much spatial structure at 
all wavelengths. Consequently, the differences in 
observed line profiles almost certainly arise from 
differences in the spatial region sampled due to 
slit sizes and angles as well as seeing conditions. 
Similarly, Winge et al. observed a [Fe II] /Pa/3 flux 
ratio of ~0.7, in comparison to our nuclear ratio 
of 0.40±0.03. Since we observe a rapidly changing 
flux ratio with position, the difference between the 
two measurements of the nuclear ratio may also be 
easily ascribed to the different spatial region sam- 
pled. 

Filippenko & Sargent (1985) report a faint 
broad (Full Width at "Zero Intensity" of 6000 



19 



km s _1 ) Ha line on the nucleus, while Shields & 
Filippenko (1988; 1996) report the detection of 
a broad line in scattered Ha flux within 5" of the 
nucleus. Figures 13 and 14 show no convincing ev- 
idence of a broad Pa/3 or Br7 line. The upper limit 
on the flux of a broad (FWHM- 2700km s" 1 ) Pa/3 
line on the nucleus is < 3 x 10~ 15 erg cm~ 2 s -1 , 
corresponding to an upper limit in its equivalent 
width of 1.6 x 10 _3 /im. 

Petitjean & Durret (1993) report on optical 
spectroscopy in slits oriented at a number of po- 
sition angles, including two angles, 0° and 60°, 
which bracket our 30° slit. To the north and east, 
within 5" of the nucleus, the optical lines show a 
shallow gradient to the red with all line velocities 
within 100 km s _1 of the systemic velocity, con- 
sistent with the infrared emission lines. Optically, 
this gradient continues to the south and west, with 
lines more distant from the nucleus at lower veloc- 
ity, in contrast to the increasing velocity seen in 
infrared Br7 and [Fe II] emission. 

Based on the offset between the centroid of the 
IR continuum image of NGC 4388 obtained with 
the spectrometer in imaging mode and a star from 
the Hubble Guide Star Catalog 3' away, the peak 
of the infrared continuum emission corresponds to 
the brighter, northeastern radio peak of two peaks 
separated by ~ 2" in the 4.86 GHz map of Hummel 
& Saikia (1991). As with Mk 1066 (section 3.1), 
there is a spatial correspondence in NGC 4388 be- 
tween enhanced [Fe II] emission and extended ra- 
dio emission. To the southwest close to the second 
peak of radio emission (Hummel & Saikia 1991), 
the [Fe II] /Pa/3 ratio increases to values > 1.5 (see 
Figure 15), whereas within 1" of the nucleus the 
[Fe II] /Pa/3 ratio is < 0.5. To the northeast, the 
[Fe II] /Pa/3 ratio is also higher than it is on the 
nucleus, but is still not as large as the values to 
the southwest. Additionally, to the southwest the 
[Fe II] line profiles are most disparate from the 
Pa/3 lines profiles, where the peak of the [Fe II] 
emission is redshifted relative to the peak of the 
Pa/3 emission. 

3.3.4- [Fe II] an d Hi emission mechanisms 

The velocity field of NGC 4388 is known to be 
complicated. There is a rotating disk with ma- 
jor axis at position angle <~ 60° and inclination 

70° (Petitjean & Durret 1993). However, in 

addition to this, there is blueshifted gas out of the 



plane of the galaxy to the north and south, as well 
as other localized sites within 20" where the veloc- 
ities of the gas are discrepant from pure rotation 
(Corbin, Baldwin, & Wilson 1988). 

In the infrared spectra, we are probably observ- 
ing at least two kinematic components within 4" of 
the nucleus of NGC 4388. One kinematic compo- 
nent is the rotating disk, which is responsible for 
the velocities of the lines northeast of the nucleus, 
and for the velocity of the H2 line southwest of the 
nucleus (Figure 17). These velocities match most 
closely the velocities of the optical lines which Pe- 
titjean & Durret (1993) identify with the rota- 
tion curve of the galaxy. In addition to this rotat- 
ing component, there may be a receding outflow 
to the southwest, associated with the radio emis- 
sion (Hummel & Saikia 1991). The shoulders 
associated with this outflow pull the Pa/3, Br7, 
and [Fe II] line centroids to higher velocities than 
would be the case from pure rotation. The veloc- 
ity difference is greater for [Fe II] than it is for 
Pa/3 because of the higher [Fe II] /Pa/3 ratio in the 
outflow than in the disk. This receding outflow 
has a higher [Fe II] /Pa/3 ratio than the gas in the 
disk of the galaxy. Colina et al. (1987) observed 
red components in the [OIII] line profiles along 
a 23° slit; to the southwest, the mean [OIII] ve- 
locities are redshifted from the peak [OIII] veloc- 
ities by 50-80 km s _1 . The Br7 line southwest of 
the nucleus probably includes contributions from 
both the disk and outflow components. The ve- 
locities of Br7 here are redshifted relative to Pa/3, 
and relative to optical emission lines that corre- 
spond to the rotation of the disk. Additionally, 
the Br7/Pa/3 ratio is higher to the southwest than 
it is to the northeast (although it is highest on the 
nucleus). This suggests that the emission associ- 
ated with the outflow is viewed through substan- 
tial extinction in the disk of the galaxy. 

The [Fe II] emission also increases off the nu- 
cleus along the 120° slit which is perpendicular 
to the putative outflow axis. This increased emis- 
sion does not correspond to notable differences in 
emission at other wavelengths. It is likely that 
this strong [Fe II] emission is associated with the 
disk of the galaxy, as it shares a similar (although 
not identical) velocity gradient with the other lines 
(Figure 18). 

The H 2 emission appears to be originating pri- 
marily in the disk of the galaxy, as its velocity as 
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plotted in Figures 17 and 18 follows closely the 
Pa/3 and optical line velocities identified with the 
rotating disk (Petitjean & Durret 1993). The 
H2/Br7 ratio diminishes away from the nucleus 
along the direction of the radio emission, suggest- 
ing that the H 2 emission is not associated with 
the radio jets or the outflow. Similar to Mk 1066, 
the H 2 /Br7 ratio shows a striking increase along 
the slit perpendicular to the outflow axis. Along 
this direction, the H2 emission appears to share 
a similar spatial structure to the [Fe II] emission. 
The ratios of both lines to hydrogen recombina- 
tion lines (Figure 15) and the velocity curve of 
both lines (Figure 18) along a position angle of 
120° share the same shape (although the H 2 emis- 
sion increases faster relative to Br7 than the [Fe II] 
emission does relative to Pa/3) . The mechanism of 
excitation of the [Fe II] and H 2 lines in the disk of 
NGC 4388 at radii of >150pc are probably related. 

3.3.5. Summary 

1. Strong Pa/3, Br7, [Fe II], and H 2 emission 
is seen from NGC 4388. These lines are ob- 
served to be extended by ~ 8" along the 
direction of greatest extent of optical [OIII] 
line emission. Infrared line emission is seen 
extended by ~ 5" perpendicular to this di- 
rection. 

2. We identify a line at the rest wavelength 
of 1.252/im with the infrared coronal line 
[S IX]. This line is only marginally spatially 
resolved at a position angle of 30°. On the 
nucleus, its flux is 1/6 the flux of Pa/3. 

3. The [Fe II] /Pa/3 flux ratio is lowest on the 
nucleus with a value of 0.4. It rises to val- 
ues > 2 southwest of the nucleus where the 
velocity centers of the Pa/3 and [Fe II] lines 
are most different, spatially coincident with 
a secondary peak of radio emission. The 
H 2 /Br7 flux ratio is highest (>2-3) away 
from the nucleus along the 120° slit, per- 
pendicular to the axis of the extended [OIII] 
emission. Along this direction, the H 2 /Br7 
and [Fe II] /Pa/3 ratios both increase away 
from the nucleus. Along the 30° slit, the two 
line ratios show opposing trends, with the 
H 2 /Br7 ratio increasing and the [Fe II] /Pa/3 
ratio decreasing with distance from the nu- 
cleus. 



4. Northeast of the nucleus, the infrared lines 
in NGC 4388 show structure which is simi- 
lar to the optical line measured by Petitjean 
& Durret (1993). Southwest of the nucleus, 
the Br7 and [Fe II] lines are redshifted rela- 
tive to Pa/3 by 50-150 km s" 1 . 

5. In some cases where different lines have dif- 
ferent velocity centroids, there is substruc- 
ture suggestive that the same velocity com- 
ponents may be present in every line with 
different strengths. Specifically, southwest 
of the nucleus where the [Fe II] emission is 
redshifted relative to the Pa/3 emission, there 
are enhanced blue wings in [Fe II] and en- 
hanced red wings in Pa/3. 

6. A possible explanation for the velocity struc- 
tures seen in the infrared lines is that there 
are at least two kinematical components. 
Most velocities are governed primarily by a 
rotating disk. There may additionally be 
an obscured redshifted outflow southwest of 
the nucleus. The coincidence of strong, red- 
shifted [Fe II] emission to the southwest may 
be evidence of shock-excited [Fe II] emitting 
gas associated with this outflow. 

3.4. Mk 3 

Mk 3 is an SO galaxy at a systemic velocity 
of 4046 km s" 1 (Whittle et al. 1988). This gives 
a calculated distance of 54 Mpc, which in turn 
yields a spatial scale of 262 pc arcsec -1 . Nar- 
rowband optical imaging with the HST in several 
lines include [OIII] and Ha shows that the NLR of 
the galaxy is S-shaped (Capetti et al. 1995). On 
small spatial scales of <~ 1", the NLR is at a po- 
sition angle of 86°. On larger spatial scales of a 
few arcseconds, the apparent position angle of the 
NLR is 113° (Pogge & de Robertis 1993). Ra- 
dio observations at 14.9 GHz show a linear struc- 
ture oriented at a position angle of 86° (Kukula 
et al. 1993). There are two radio peaks, the 
stronger one 1.1" west of the nucleus, the weaker 
one 0.4" east of the nucleus. Higher resolution 
maps of the source at 5 GHz resolve the linear 
structure into a number of components. In partic- 
ular, the eastern lobe is resolved into three compo- 
nents aligned along the axis of the source, together 
with a third component coincident with the core 
of the galaxy. Mk 3 is also quite bright in X-rays, 
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with an X-ray luminosity of 3 x 10 42 erg s 1 



(with 



Ho = 75 km s -1 Mpc -1 ), and a very high absorb- 
ing column of 6 x 10 23 cm~ 2 (Awaki et al. 1990). 
The 0.1-2 keV X-ray source is unresolved with the 
ROSAT HRI on a scale of - 4" (Morse et al. 1995). 

The reduced 1.29^m and 2.21/xm spectra of 
Mk 3 along a 113° slit are plotted in Figure 20, 
in spatial bins along the slit separated by 1". The 
position angle of the 113° slit was chosen to align 
with the major axis of the [OIII] distribution on 
scales of - 1" (Mulchaey et al. 1996). The [Fe II] 
and Pa/3 emission lines are visible in the 1.2/xm 
spectra, and the H2 and Br7 emission lines are 
visible in the 2.2^m spectra. Note that for these 
observations, the Br7 line was centered in the K- 
band spectrum so that the H2 line is near the blue 
edge of the spectrum. This may lead to inaccura- 
cies in placing the continuum on the blue side of 
the H 2 , resulting in underestimates of the flux of 
the H 2 line and any blue wings in the profile of the 
H 2 line. 

The [Fe II] , Pa/3, H 2 , and Br7 lines are extended 
along the slit by < 7" (1.8 kpc). The equivalent 
widths of the lines on the nucleus, given in Table 5 
in Appendix A, are larger by a factor of > 2 than 
those seen in most Seyfert 2 nuclei (Knop et al. 
2000). 

Most of the extended line emission is charac- 
teristically narrow (FWHM< 500 km s _1 ), yet on 
the nucleus the infrared lines are quite broad (see 
Figure 21). The FWHM of the nuclear lines are 
Br 7 -700 km s"\ Pa/3 -1000 km s"\ and [Fe II] 

— 1150 km s _1 . These lines are also markedly 
asymmetric, with strong blue wings. The peak 
of the Pa/3 line on the nucleus is at a velocity 

— 200 km s _1 greater than the centroid of the 
full line provile. Although [S IX] is not obvious, it 
could easily be hidden within the broad blue wings 
of the [Fe II] line profile 

While the observed width of the Pa/3 line on the 
nucleus is comparable to the FWHM of broad Pa/3 
lines seen in some Seyfert 1 galaxies (e.g., Knop et 
al. 2000), the observed line profiles in Mk 3 prob- 
ably do not indicate the detection of a canonical 
broad line region, since [Fe II] and Pa/3 lines are 
both broad and have similar profiles. Because the 
critical density of [Fe II] is only n c — 10 5 cm~ 3 , 
the broad [Fe II] emission cannot be produced in 
a high density (n c — 10 8 — 10 10 cm~ 3 ) broad line 
region. Some or all of the high velocity Pa/3 line 
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Fig. 20.— Spectra of Mk 3 along the 113° slit. 
Each spectrum is within a rectangular beam which 
is 1.0" along the slit by the width of the slit (— 
0.6"). Spatial bins are adjacent. The distance of 
each spatial bin from the nucleus in an easterly 
direction is indicated in the spectrum. Flux units 
(Fa) arc arbitrary. 
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Fig. 21. — Profiles of nuclear lines in Mk 3. The 
lines are normalized to the same peak flux. The 
solid line is Pa/3, the dotted line [Fe II], and the 
dashed line Br7. 



flux is therefore probably coming from the NLR. 

3.4.1. Flux Ratios 

Table 5 in Appendix A lists the equivalent 
widths and flux ratios of the lines observed 
in Mk 3. Figure 22 shows the flux ratios of 
[Fe II] /Pa/3 and H 2 /Br7 as a function of position 
along the slit. The [Fe II] /Pa/? ratio changes by 
a factor of nearly 4 along the slit. This change is 
a monotonic increase from east to west, although 
it is poorly determined 3" west of the nucleus due 
to the low signal to noise ratio in the Pa/3 line. 
The H2/Br7 flux ratio, close to 0.3 — 0.4 in all but 
one spatial bin, is relatively low for Seyfert galax- 
ies. The H2/Br7 flux ratio is anomalously large 2" 
west of the nucleus. The strong, double-peaked H2 
emission responsible for this enhanced flux ratio 
is visible in Figure 20. The higher velocity peak is 
unresolved, and is at a velocity of 450 km s^ 1 rel- 
ative to the systemic velocity of Mk 3. This peak 
is not well detected in every individual image that 
goes into the summed spectrum of Mk 3, and so 
we suspect that this peak is an artifact associated 
with poor background subtraction. Because H 2 at 



the redshift of Mk 3 falls right at the wavelength of 
an atmospheric OH line, the spatial bins in which 
the detection of H2 is of low significance may be 
containmated by poor subtraction of this line. If 
this is the case, then the H 2 /Br7 ratio would be 
< 1.5, which is still higher than in the other spa- 
tial bins, although the caveat about background 
subtraction applies to the full H 2 profile in this 
spatial bin. 
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Fig. 22. — Mk 3 flux ratios as a function of posi- 
tion along the 113° slit. Positive positions are to 
the east. Fluxes are integrated over the entire line 
profile in each spatial bin. 



3.4-2. Velocity Structure 

Figure 23 shows position versus velocity plots 
for each of the four emission lines seen in Mk 3. 
A smoothed continuum has been subtracted from 
the data in each case. Table 9 in Appendix A lists 
the velocity centroids and total full widths at half 
maximum of the lines in each of the 1" spatial bins 
for which spectra are plotted in Figure 20. 

To the east, the Pa/3 line is narrow and stronger 
than any narrow Pa/3 emission to the west. A very 
shallow velocity gradient at the limit of our reso- 
lution is marginally visible in this narrow compo- 
nent. The Pa/3 line is redshifted by <~ 100 km s _1 
> 3" east of the nucleus. The Br7 line is not as 
well detected as the Pa/3 line, but appears to share 
a similar velocity structure, although the shallow 
velocity gradient is not apparent in Br7. The Br7 
line is slightly blueshiftcd relative to the Pa/3 line 
east of the nucleus. The Pa/3 and Br7 line peaks 
appear to share a shallow velocity decrease from 
east to west. The total blueshift of the lines is 
< 100 km s _1 west of the nucleus. 

The velocity structure of the [Fe II] line to the 
east is different from that of the hydrogen recombi- 
nation lines. To the east of the nucleus and on the 
nucleus, the peak of the [Fe II] line is blueshiftcd 
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Fig. 23. — Position versus velocity plots for Mk 3. 
Position is along a 113° slit. Velocity is relative 
to 4046 km s _1 , the systemic velocity of Mk 3 
(Whittle et al. 1988). Contours are logarithmic, 
separated by a factor of 



by ~ 100-200 km s" 1 relative to the Pa/3 and Br 7 
lines. West of the nucleus, the velocity of the 
[Fe II] line is closer to that of the hydrogen re- 
combination lines. The H2 line is not confidently 
detected east of the nucleus. West of nucleus by 
2", the H 2 line shows the aforementioned peculiar 
double-peaked shape (see Figure 20). 

On the nucleus, the [Fe II], Pa/3, and Br 7 lines 
all have qualitatively similar profiles. All three 
lines show strong enhanced emission to the blue of 
the peak. However, as is visible in Figure 21, the 
profiles differ in detail. In Br 7 , the blue wing is not 
as strong relative to the peak as it is for Pa/3 and 
[Fe II] . The peak of the [Fe II] line profile is shifted 
to the blue of the Pa/3 line profile by 80 km s _1 , 
and the FWHM of [Fe II] is about 100 km s" 1 
larger than the FWHM of Pa/3, apperently due to 
enhanced emission on the blue wing of [Fe II] (see 
Figure 21). 

While the profiles to the east of the nucleus are 
narrow, Figure 23 clearly shows that 1-2" west of 
the nucleus the lines are broader, with the FWHM 



of Pa/3 line equal to 500-600 km s _1 . Figure 24 
shows the profiles of the Br 7 , Pa/3, and [Fe II] lines 
1" west of the nucleus. Two features are apparent 
from this plot. First, the red wing of the Br 7 
line is stronger relative to the line peak than the 
red wing of Pa/3. This may indicate that there is 
receeding gas obscured behind dust, as the longer 
wavelength Br 7 is less affected by extinction than 
Pa/3. Second, the [Fe II] /Pa/3 ratio is largest with 
a value of <~ 2.2 at a redshifted velocity near <~ 
150 km s _1 . The blue hump on the Pa/3 and [Fe II] 
line profiles has an [Fe II] /Pa/3 ratio of about 1, 
close to the nuclear flux ratio (see section 3.4.1). 

3.4-3. Comparison with Other Observations 

The velocity gradient seen in the off-nuclear 
narrow Pa/3 line is generally consistent with the 
rotation curve seen in optical H/3 and [OIII] emis- 
sion by Whittle et al. (1988). However, the optical 
slit position angle is 86°, while the infrared data 
presented here was obtained through a slit at a po- 
sition angle of 113°. East of the nucleus, the [OIII] 
line peaks at 120 km s _1 , while west of the nucleus 
it peaks at -120 km s _1 , similar to the limiting off- 
nuclear velocities of ±100 km s _1 in the infrared 
lines. Whittle et al. (1988) also saw an increase 
in the FWHM of the optical H/3 and [OIII] lines 
near the nucleus, similar to that seen in the Pa/3, 
Br 7 , and [Fe II] lines. This broader emission is 
blueshifted relative to the peak of the lines, again 
consistent with that seen in the infrared lines. In 
addition, a red wing in the [OIII] and H/3 profiles 
is similar to the enhanced red emission seen in 
the Br 7 and [Fe II] lines in Figure 24, which was 
also strongest 1" west of the nucleus. A Gaus- 
sian decomposition of the optical lines placed the 
relative velocity of the redshifted component at 
~ 300 km s ; our Gaussian decomposition of 
the infrared lines places the velocity of the red- 
dened component of the infrared lines at some- 
what higher velocity, closer to ~ 500 km s _1 , al- 
though this is uncertain given the complex blended 
line profiles (see Figure 24). 

The spatial location of this red wing is coinci- 
dent with the stronger of two radio peaks seen at 
2cm (Kukula et al. 1993). The relatively high ex- 
tinction seen for this reddened component in the 
infrared Br 7 and Pa/3 lines is consistent with the 
picture of this radio peak as a lobe coincidental 
with a receding outflow of ionized gas which is ob- 
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Fig. 24. — Line profiles of Br7, [Fe II], and Pa/?, 
and the log of the [Fe II] /Pa/3 flux ratio as a func- 
tion of velocity, 1" west of the nucleus of Mk 3 at 
a position angle of 293°. Points with uncertainty 
greater than 20% have been omitted from the lat- 
ter plot. 



scured behind the plane of the galaxy. The weaker 
radio lobe is 0.4" east of the nucleus along the 
same position angle. Whittle et al. (1988) do not 
identify any component in the optical spectra with 
the eastern radio lobe. However, the enhanced 
blue wings seen on the nucleus in both the in- 
frared and [OIII] emission lines may be part of an 
approaching outflow associated with the eastern 
radio lobe. 

3.4-4- [F e H] an d H2 Emission Processes 

As with Mk 1066 and NGC 4388, and to a lesser 
extent NGC 2110, the [Fe II] emission in Mk 3 
shows spatial correlations with the radio emission 
which suggests that processes involving the radio 
jet's interaction with circumnuclear gas are signif- 
icantly enhancing the [Fe II] emission. West of the 
nucleus at the position of the western radio peak, 
the [Fe II] line shows enhanced redshifted emis- 
sion, which is also seen in the Pa/? and Br7 lines 
as well as in the optical lines (Whittle et al. 1988). 
At the position of this radio peak west of the nu- 
cleus, the [Fe II] /Pa/3 ratio is higher than it is to 
the east of the nucleus where the radio emission is 
much weaker. This enhanced redshifted emission 
1" west of the nucleus also has a larger Br7/Pa/3 
flux ratio, suggesting higher extinction to the re- 
gion which is the source of the enhanced [Fe II] 
emission. This is consistent with the picture of 
the region as a receding outflow which is partially 
obscured by the galactic disk. 

The eastern radio peak is too close (0.4") to 
the nucleus to distinguish infrared line emission 
spatially coincident with it from emission directly 
from the nucleus. However, if the two knots of ra- 
dio emission represent the radio lobes associated 
with a biconical outflow, the enhanced blue wings 
in the infrared and optical lines on the nucleus may 
in fact be correlated with the second radio lobe. 
Higher resolution 5 GHz maps of Mk 3 (Kukula 
et al. 1993) resolve the eastern lobe into a number 
of individual linearly spaced components. The in- 
frared lines near the nucleus show a complicated, 
evidently multi-component structure in velocity 
(see Figure 20). The multi-component spatial na- 
ture of the radio emission and the multiple velocity 
components in the infrared line emission suggest 
the identification of the broad blue wings on the 
infrared lines near the nucleus with the eastern 
radio lobe. The [Fe II] emission on the nucleus 
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is blueshifted relative to the Pa/3 emission, sug- 
gesting that the [Fe II] /Pa/3 flux ratio may be en- 
hanced at the position of the eastern radio peaks. 
Given this spatial correspondence, it is likely that 
the [Fe II] emission is enhanced by fast shocks re- 
sulting from the interaction of the radio jet and 
the circumnuclear gas. 

Except for the anomalously high value 2" west 
of the nucleus, the H 2 /Br7 ratio is lower in Mk 3 
than it is in most Seyfert galaxies. West of the nu- 
cleus, the narrow Pa/3 and perhaps Br7 emission 
may be associated with the rotation of a galac- 
tic disk, based on similarities with optical lines 
(Whittle et al. 1988). At these positions, the H 2 
emission is suppressed relative to the Br7 emis- 
sion, indicating that Mk 3 continues relatively lit- 
tle warm molecular gas. 

3.4-5. Summary 

1. At a position angle of 113°, the infrared line 
emission is observed to be extended by < 7" . 

2. There is a FWHM~1000 km s" 1 blue wing 
visible on the nucleus of Mk 3, most most 
prominent in [Fe II] and Pa/3 but also seen 
in Br7. There is a red wing on the line pro- 
files, arising from a region of relatively high 
extinction, <~ 1" west of the nucleus. This 
red wing corresponds to a rcdshifted compo- 
nent seen in optical emission (Whittle et al. 
1988). Both wings are spatially consistent 
with the position of two radio peaks (Kukula 
et al. 1993). 

3. The Pa/3 and Br7 narrow line peaks may fol- 
low a barely resolved smooth velocity gradi- 
ent which lies at +100 km s _1 cast of the 
nucleus and —100 km s -1 west of the nu- 
cleus. East of the nucleus, the Br7 line is 
marginally blueshifted relative to Pa/3. The 
[Fe II] emission is blueshifted relative to the 
Pa/3 emission both east of the nucleus, and 
significantly so on the nucleus. 

4. The [Fe II] /Pa/3 flux ratio monotonically in- 
creases from <~ 0.5 to ~ 2 over 5" from east 
to west across the nucleus. The H 2 /Br7 flux 
ratio is close to 0.3 — 0.4 everywhere except 
where it is anamalously high, possibly as a 
result of an artifact, 2" west of the nucleus. 



We suggest that emission associated with an 
outflow of ionized gas contributes significantly to 
the infrared line emission in Mk 3. The [Fe II] /Pa/3 
ratio is enhanced in components identified with 
this outflow and coincident with the 2cm radio 
lobes, suggesting a possible shock contribution to 
the [Fe II] emission. The low H 2 /Br7 ratio indi- 
cates that in contrast to other Sefyert galaxies, the 
circumnuclear disk is relatively depleted in warm 
molecular gas. 

4. General Discussion 

4.1. [Fe II] Emission Processes 

One obvious conclusion from the data presented 
here is that the [Fe II] /Pa/3 ratio in Seyfert galax- 
ies varies significantly with position on scales 
of < 1". Of the four galaxies observed, each 
showed a different dependence of [Fe II] /Pa/3 on 
position, whether it is monotonically increasing 
across the nucleus (Mk 3), highest at the nucleus 
(NGC 2110), lowest at the nucleus (NGC 4388), 
or close to constant with point-to-point variations 
(Mk 1066). These data show that flux ratios inte- 
grated over large nuclear beams may be of limited 
value for drawing conclusions about the excita- 
tion processes of the infrared lines. The flux ra- 
tios are observed to vary significantly on scales of 
< 10 2 pc, and to show significant velocity struc- 
ture. 

Circumnuclear starbursts and processes asso- 
ciated with the central AGN directly have both 
been implicated in the generation of [Fe II] emis- 
sion lines seen in Seyfert galaxies (e.g., Moorwood 
& Oliva 1988; Mouri et al. 1993). Our data in 
general support the active nucleus as the power 
source of [Fe II] in the Seyfert galaxies. Despite 
the differences in the details of the velocity struc- 
ture of the infrared lines, the velocity gradients 
of [Fe II] tend to match those of Pa/3 and Br7. 
The [Fe II] line can generally be kinematically as- 
sociated with other spectral features linked with 
the nuclear activity (e.g., [OIII]), or sites of en- 
hanced radio emission. The one exception is the 
extended [Fe II] emission along position angle 120° 
in NGC 4388. 

There are two primary proposed mechanisms 
for [Fe II] emission associated with the AGN it- 
self. Because the ionization potential of Fe + is 
only 16.2 eV, the [Fe II] line will be strong in par- 



26 



tially ionized regions. If clouds illuminated by the 
Seyfert nucleus have a high enough column den- 
sity, only the side of the cloud facing the nuclear 
ionizing source will be fully ionized. Because the 
ionization cross section for high-energy X-rays is 
relatively small, X-rays can penetrate deeply into 
these dense clouds, creating large warm partially 
ionized regions suitable for [Fe II] emission (Gra- 
ham et al. 1990). Fast shocks can also produce 
partially ionized regions heated to 6000 K (Shull & 
Draine 1987). Models indicate that these shocks 
may also destroy dust grains (Seab & Shull 1983). 
Because much of the iron in the local interstel- 
lar medium is locked up in dust grains, the de- 
struction of those grains would enhance the gas 
phase Fe abundance, thereby increasing the in- 
frared [Fe II] emission. 

Simpson et al. (1996b) report a slightly stronger 
correlation between infrared [Fe II] and 6cm ra- 
dio emission than between Pa/3 and 6cm emission. 
Based on the difference in the correlations, Simp- 
son et al. conclude that <~ 20% of the [Fe II] emis- 
sion in Seyfert galaxies is due to the action of fast 
shocks, and that the bulk of it is due to X-ray 
photoionization. 

The data presented here for Mk 1066, 
NGC 2110, NGC 4388, and Mk 3, combined with 
the calculations of Simpson et al, provide spa- 
tial and kinematic evidence that indicate that 
X-ray photoionization is responsible for most of 
the [Fe II] emission, but that the [Fe II] emission 
is probably enhanced in some places by a contri- 
bution from shocks. In most cases the [Fe II] emis- 
sion tends to be spatially associated with the Pa/3 
emission, and the primary component of [Fe II] 
often shares the same velocity structure as the 
Pa/3 line. Thus, most of the [Fe II] emission prob- 
ably comes from gas physically and kincmatically 
associated with the NLR, as traced by the Pa/3 
emission. Except for the nucleus of NGC 2110, 
the flux ratios observed at various positions in the 
galaxies are within the range of those seen in the 
spectra of Simpson et al. 

While the bulk of the [Fe II] appears to be pho- 
toionized, there are clear signatures of local en- 
hancements of the [Fe II] emission from shocks as- 
sociated with the expanding radio plasma. In all 
four of the galaxies, spectral features were identi- 
fied in the [Fe II] emission line profile which oc- 
cured at or near the position of peaks in the non- 



thermal radio emission. In many cases, these fea- 
tures showed relatively high [Fe II] /Pa/3 ratios. In 
Mk 3 and NGC 4388, the [Fe II] emission was en- 
hanced relative to the Pa/3 emission, or showed 
the largest kinematical differences from the Pa/3 
emission, at positions spatially coincident with ra- 
dio peaks. This suggests that processes associated 
with radio jets are contributing measurably to the 
[Fe II] emission. The most probable explanation 
of this is fast shocks resulting from the interaction 
of outflowing radio plasma with the circumnuclear 
gas. 

4.2. H 2 Emission Processes 

Molecular hydrogen emission can result from 
ultraviolet fluorescence of hydrogen molecules in 
photodissociation regions illuminated by pho- 
tons from OB stars or by a nuclear UV con- 
tinuum source. Alternatively, the H 2 emission 
may be excited thermally, through either slow 
(v< 25 km s" 1 ) shocks (Shull & Draine 1987), or 
through X-ray heating of dense molecular clouds 
(Draine & Woods 1990). Observations of the flux 
ratio between H 2 (u=l-0 S(l)) (A=2.1212^m) and 
H 2 (z/=2-l S(l)) (2.2471^m) rule out ultraviolet 
fluorescence as a significant mechanism in most 
Seyfert galaxies (Moorwood & Oliva 1990). 

In all four of the galaxies discussed here, we 
have identified features in [Fe II] and hydrogen re- 
combination line profiles which we associate with 
outflows or other interactions with radio lobes and 
jets. Except for NGC 2110, these features are 
weaker or invisible in the H 2 line profile; the kine- 
matics of the H 2 arc usually better explained by 
galactic rotation. In addition, in those two galax- 
ies where we have data perpendicular to the out- 
flow axis (Mk 1066 and NGC 4388), and possibly 
as well in NGC 2110 (S99), the H 2 is significantly 
more extended than the atomic hydrogen recombi- 
nation lines, with the H 2 /Br7 flux ratio increasing 
by a factor of ^4 or more at 200-300 pc away from 
the nucleus. 

Although X-ray heating of molecular clouds 
that cross into the ionization cones can give rise 
to the extended H 2 emission along the cone axes, 
it is unlikely that the nucleus is responsible for 
the enhanced H 2 emission along the axis perpen- 
dicular to the cones. According to the standard 
model of Seyfert galaxies, the radiation from the 
central source is obscured except along the axis of 



27 



the ionization cones; thus, nuclear X-rays cannot 
penetrate significantly into the gas perpendicular 
to the cone axis to produce the extended H2 emis- 
sion. While these X-rays would heat the inner face 
of the obscuring molecular torus, the size of that 
torus is believed to be on the order of a few par- 
sees, so this cannot explain the extended H2 emis- 
sion that we observe. Also, if molecular hydrogen 
emission does emerge from photodissociation re- 
gions heated by X-rays from a central source, we 
would expect strong Br7 emission illuminated by 
ultraviolet radiation from the same central source. 

Shocks driven directly or indirectly by the nu- 
cleus are another possible explanation for the ex- 
citation of the H 2 emission. Slow shocks mov- 
ing through high density molecular clouds in the 
disk of the galaxy could excite molecular hydro- 
gen emission. Where the shocks encounter lower 
density clouds, they should be faster, and produce 
substantial [Fe II] emission. The [Fe II] /Pa/3 ratio 
does rise with distance from the nucleus, at least 
initially, perpendicular to the primary direction 
of extended high-excitation gas in both Mk 1066 
and NGC 4388 . Alternatively, secondary shocks 
in molecular clouds circumscribing radio outflows 
might give rise to enhanced H2 emission. This 
process might produce the differing [Fe II] and H 2 
spatial morphologies observed in NGC 1068 by 
Blietz et al. (1994). One problem with this sce- 
nario is that, in contrast to the axis of ioniziation 
cones where outflowing radio plasma may drive 
shocks into the surrounding gas, there is no ob- 
vious mechansim to drive shocks to distances of 
~1 kpc perpendicular to the cone axis. 

The relatively large extent of the H2 emission 
perpendicular to the ionization cones in Mk 1066 
and NGC 4388 suggest that there are local source 
of heating which is giving rise to much of the H2 
emission. In both galaxies, there may be contri- 
bution to the H2 emission from young stars heat- 
ing molecular clouds. While this may contribute 
to the H2 emission, it is unlikely to be the only 
source, as the H 2 /Br7 ratios observed in the ex- 
tended H2 gas tend to be higher than the val- 
ues < 0.5 typically observed in starburst galaxies 
(Moorwood & Oliva 1988, 1990; Kawara et al. 
1988). NGC 4388 has a bar oriented at position 
angle - 90° (Corbin, Baldwin, & Wilson 1988). 
Dynamical processes associated with this bar may 
be creating shocks which give rise to the observed 



H2 emission. In Mk 1066, it is also possible that 
the bulk motions of clouds which do not move com- 
pletely with the rotation of the disk are colliding 
with other clouds, creating slow shocks which ex- 
cite the cold molecular gas. In Mk 3, the lower 
H 2 /Br7 ratio combined with the narrower H 2 line 
profile suggest that local star formation activity 
may be able to explain most of the H 2 emission, 
with no need for a contribution from outflowing 
gas associated with the nucleus. By contrast, in 
NGC 2110, where wings associated with the out- 
flow are seen in H 2 , nuclear X-ray heating may be 
a more significant contributor. 

5. Summary and Conclusions 

1. We have observed strong infrared Pa/3, Br7, 
[Fe II], and H2 line emission which is spa- 
tially extended on scales of a few hundred 
parsecs in Mk 1066, NGC 2110, NGC 4388, 
and Mk 3. 

2. We have detected the near infrared coronal 
line [S IX] (A=1.2525/im) in the nucleus of 
NGC 4388. In contrast to the other infrared 
lines, this line is only marginally spatially ex- 
tended. There is no apparent [S IX] emission 
from Mk 1066, NGC 2110, or Mk 3. This is 
similar to what is seen optically, that high 
ionization species are visible in some but not 
all Seyfert galaxies. 

3. While often showing similar velocity trends 
over the central few arcseconds, the [Fe II], 
H2, and H lines have profiles which dif- 
fer in detail from each other over the full 
spatial range probed by these observations. 
The peaks of lines are shifted relative to 
each other at some spatial positions in some 
galaxies, and some species show enhanced 
blue and red wings relative to other species. 
These differences are usually due to the ef- 
fects of extinction, an enhancement to the 
[Fe II] emission associated with other indica- 
tors of the action of an outflow, or a possible 
local source of heating of H 2 . 

4. The [Fe II] /Pa/3 and H 2 /Br7 flux ratios can 
vary dramatically on scales of 10 2 pc near 
the nuclei of Seyfert 2 galaxies. The details 
of the spatial structure of the [Fe II] /Pa/3 ra- 
tio frequently support the notion that par- 
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tially ionized regions created by X-ray pho- 
toionization from the nucleus, or shocks from 
the interaction of outflowing gas with ambi- 
ent gas, or both, are enhancing the circum- 
nuclear [Fe II] emission. The spatial varia- 
tion of the H 2 /Br7 ratio is again frequently 
ascribable to local heating of H 2 . 



5. The [Fe II] emission appears to be associ- 
ated with the nuclear activity, not requiring 
the invocation of a circumnuclear starburst 
in any of these four galaxies. The [Fe II] can 
be kinematically and spatially linked to pro- 
cesses dominated by the active nucleus. In 
all of these galaxies, there are correlations 
between features in the [Fe II] emission and 
both radio jets and optically identified out- 
flows, suggesting that some [Fe II] emission 
is produced by fast shocks resulting from in- 
teractions between radio jets and circumnu- 
clear gas. 

6. The H 2 emission, although peaked on the 
nucleus, is relatively strong compared to hy- 
drogen recombination and [Fe II] emission 
perpendicular to the axis of ionization cones 
or the major axis of high-excitation gas in 
Mk 1066 and NGC 4388. The spatial distri- 
bution and line profile of H 2 generally differs 
from that of [Fe II]. A local source such as 
bar-driven shocks or circumnuclear star for- 
mation may be partially responsible for the 
H 2 emission in these galaxies. 

The authors would like to thank the staff and 
night assistants at Palomar Observatory for as- 
sistance in completing these observations. Addi- 
tionaly, they would like to thank James Larkin for 
discussion and assistance in the earlier phases of 
the project. 
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A. Measurements of Equivalent Widths, Fluxes, and Line Ratios 

In this appendix are tables with data from the infrared spectra of Mk 1066, NGC 2110, NGC 4388, and 
Mk 3. Tables 2 through 5 list the equivalent widths for the integrated line profile for each species observed 
in each spatial bin. They additionally list relevant flux ratios. Tables 6 through 9 give the centroids and line 
widths (FWHM) for each species in each spatial bin. 

Table 2: Mk 1066 Equivalent Widths and Flux Ratios 







Equivalent Width (10~Vm) 


Flux Ratios 


Pos a 














(") 


[Fc II] 


Pa/3 


H 2 


Br7 


[Foil] /Pa/3 


H 2 /Br 7 








Mk 1066, 135° slit 






2.67 


6.7±4.5 


17.1±12.4 






0.43±0.14 




2 


4.7±1.7 


5.8±1.9 


4. 0±1. 9 


<4.3 


0.79±0.21 


>0.87 


1.33 


18.4±3.2 


15.3±2.6 


9.4±2.0 


10.3±2.7 


1.17±0.08 


0.86±0.09 


0.67 


23.0±2.6 


29.0±2.6 


12.6±1.7 


18.7±2.8 


0.75±0.03 


0.63±0.05 





20.8±1.9 


27.4±2.1 


13. Oil. 2 


14.0±1.5 


0.75±0.04 


0.89±0.06 


-0.67 


39.0±3.8 


33.0±4.6 


17.8±2.6 


20.6±3.8 


1.20±0.04 


0.85±0.06 


-1.33 


21.6±3.4 


23.2±4.6 


9.0±2.7 


16.7±6.3 


0.99±0.05 


0.56±0.09 


-2 


5.3±2.2 


6.7±2.5 


<4.4 


11.5±6.1 


0.85±0.17 


<0.40 








Mk 1066, 


45° slit 






1.33 


5. Oil. 3 


4.2±1.1 


13.0±5.3 


<2.7 


1.30±0.32 


>4.41 


0.67 


11.2±1.6 


10.3±1.3 


12.6±2.7 


6.7±1.1 


1.06±0.09 


1.74±0.21 





17.2±1.4 


29.1±2.7 


11.5±1.5 


13.7±1.4 


0.57±0.03 


0.76±0.06 


-0.67 


15.9±2.6 


17.9±2.7 


12.6±2.3 


10.0±2.0 


0.87±0.07 


1.21±0.10 


-1.33 


5.2±1.4 


8.9±2.6 


9.0±4.0 


<4.7 


0.64±0.16 


>1.83 


-2 






5.6±3.7 


<5.3 




>0.91 



a Positive positions are towards the position angle of the slit as quoted. 



Table 3: NGC 2110 Equivalent Widths and Flux Ratios 







Equivalent Width (10~ 4 




Flux Ratios 


Pos a 














(") 


[Fc II] 


Pa/3 


H 2 


Br-y 


[Fc II] /Pa/3 


H 2 /Br 7 


3 


3.5±1.0 


3.3±0.8 


3.8±1.3 


2.1±0.9 


1.09±0.25 


1.89±0.78 


2 


6.2±1.3 


2.9±0.7 


7.2±1.4 


1.9±0.6 


2.19±0.42 


4.08±1.08 


1 


15.4±1.4 


3.1±0.5 


6.2±0.8 


2.0±0.4 


5.02±0.48 


3.23±0.53 





24.7±1.6 


3.1±0.4 


4.8±0.4 


1.6±0.3 


8.11±0.78 


3.17±0.63 


-1 


22.5±2.3 


3.2±0.5 


6.6±0.7 


2.0±0.4 


7.36±0.74 


3.59±0.69 


-2 


11.7±2.3 


2.8±0.6 


8.4±1.7 


1.4±0.4 


4.43±0.53 


6.50±1.50 


-3 


8.7±2.0 


3.1±0.9 


5.5±1.2 


2.8±1.0 


2.93±0.54 


2.14±0.51 


-4 


5.8±2.8 


<2.7 


3.6±1.5 


<9.7 


>2.24 


>0.39 



a Positive positions are toward the southeast at a position angle of 160°. 
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Table 4: NGC 4388 Equivalent Widths and Flux Ratios 







Equivalent Width (10~Vm) 






Flux Ratios 




Pos a 


















f"1 

I ) 


TFp TT1 


Pa/3 


H 2 


Br7 


[S IX] 


fFe TTl /Pa/3 


112 / 1J1 / 


IS TXl /Pa/3 










NGC 4388, 30° slit 








4 


22.8±19.1 


19.0±14.9 


<7.4 


13.0±12.8 




0.99±0.25 


<0.67 




3 


15.5±9.5 


18.2±9.1 


<7.0 


13.3±9.9 




0.71±0.17 


<0.57 




2 


13.4±6.0 


26.1±7.0 


6.5±3.1 


12.7±6.7 




0.47±0.08 


0.56±0.13 




1 


7.7±1.9 


25.7±5.5 


13.6±4.1 


12.7±3.2 


3.6±1.0 


0.28±0.03 


1.05±0.15 


0.13±0.02 





11.2±1.7 


26.3±4.2 


9.5±1.2 


10.6±1.5 


4.5±0.8 


0.40±0.03 


0.88±0.08 


0.16±0.02 


-1 


18.8±4.6 


32.2±6.8 


8.4±1.6 


17.1±3.9 


2.4±0.8 


0.57±0.04 


0.51±0.06 


0.07±0.01 


-2 


31.0±9.4 


26.8±7.4 


6.3±2.7 


14.7±6.5 




1.13±0.11 


0.46±0.09 




-3 


29.7±11.8 


17.8±7.1 


6.1±3.6 


13.6±10.5 




1.64±0.22 


0.49±0.14 




-4 


12.5±6.6 


4.8±2.7 


<6.9 


<11.7 




2.52±0.84 














NGC 4388, 


120° slit 








2 


9.0±2.5 


4. Oil. 3 


12.4±5.3 


<4.9 




2.18±0.53 


>2.29 




1 


12.7±2.1 


16.8±3.1 


14.1±2.9 


9.4±2.0 




0.72±0.07 


1.39±0.18 







11.5±1.4 


24.1±2.7 


8.4±1.1 


11.6±1.8 


4.4±0.6 


0.46±0.03 


0.67±0.07 


0.17±0.01 


-1 


10.1±2.0 


19.9±3.5 


19.6±4.6 


8. Oil. 7 




0.50±0.05 


2.34±0.25 




-2 


8.9±3.0 


7.2±2.2 


18.0±5.9 


3.2±1.5 




1.26±0.24 


5.23±1.60 





a Positive positions are towards the position angle of the slit as quoted. 



Table 5: Mk 3 Equivalent Widths and Flux Ratios 







Equivalent Width (10 4 




Flux Ratios 


Pos a 














(") 


[Fc II] 


Pa/3 


H 2 


Br7 


[Fell] /Pa/3 


H 2 /Br 7 


3 


11.6±4.3 


24.5±8.6 


5.8±4.2 


12.0±8.4 


0.47±0.10 


0.49±0.24 


2 


13.0±3.5 


17.7±4.7 




4. Oil. 4 


0.74±0.12 




1 


37.6±4.2 


31.3±3.5 


3.4±1.2 


11.5±2.6 


1.17±0.07 


0.29±0.11 





52.2±5.1 


39.3±3.2 


5.3±0.7 


16.7±2.1 


1.24±0.04 


0.31±0.04 


-1 


46.3±6.2 


26.4±3.0 


4.6±0.8 


12.8±2.2 


1.72±0.08 


0.35±0.05 


-2 


17.9±5.7 


9.5±2.9 


13.1±4.5 


4.4±1.6 


1.88±0.43 


2.94±0.54 


-3 


17.6±7.7 


4.3±1.7 






3.75±1.15 




a Positiv. 




3 cast towards the 


113° posit 


ion angle of tho slit. 
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Tabic 6: Mk 1066 Line Ccntroids and Widths 







[Fc IIJ 




Pa/3 


H 2 




Br-y 


Pos a 

II 


Centroid b FWHM C 
km s — ^ km s~ ^ 


Centra 
km s~ 


id FWHM 
1 kras" 1 


Centroid 
km s — 1 


FWHM 
km s — ^ 


Centroid FWHM 
km s~ ^ km s — ^ 



Mk 1066, 135° slit 



2 


.67 


-118±24 


< 154 


-233±41 


677±194 










2 




-68±25 


267±109 


-58±20 


221±83 


-72±35 


< 190 






1 


33 


-59±14 


294±53 


-77±19 


< 154 


-45±15 


328±46 


-90±19 


383±34 





.67 


-38±14 


240±32 


-67±14 


< 154 


-45±15 


321±50 


-73±18 


237±46 







-125±16 


296±41 


-20±15 


205±50 


3±13 


353±29 


-19±18 


329±27 


■0 


67 


-54±14 


289±54 


51±14 


< 154 


89±13 


223±37 


51±19 


269±32 


■ 1 


33 


-8=15 


341±62 


32±16 


167±74 


143±20 


360±78 


70±20 


308±44 


■ 2 




109±18 


218±75 


70±18 


202±65 






133±46 


715±196 










Mk 


1066, 45° 


slit 








1 


33 


-34±17 


222±73 


-50±19 


179±93 


-6±15 


491±59 









67 


-52±14 


179±52 


-20±10 


< 150 


-12±11 


392±35 


-43±17 


429±62 







-74±12 


179±38 


-41±13 


< 150 


-14±12 


192±61 


-83±17 


250±44 


■0 


67 


-54±15 


219±46 


-24±11 


< 150 


-21±10 


277±33 


-49±14 


313±51 


■ 1 


33 


-10±19 


247±83 


64±49 


< 150 


39±18 


329±82 







-2 ... ... ... ... 107±23 < 190 



a Positive position is towards the position angle of the slit as quoted. 

b Ccntroids arc relative to 3625 km s — ^. the systemic velocity of Mk 1066. 

c The resolution has been subtracted in quadrature. 



Tabic 7: NGC 2110 Line Ccntroids and Widths 







IFe IIJ 


Pa/3 


H 2 


Br^ 


Pos a 
1/ 


Centroid 6 FWHM C 
km s _1 km s _1 


Centroid FWHM 
km s -1 km s _1 


Centroid FWHM 
km s -1 km s -1 


Centroid FWHM 
taT 1 kms- 1 



3 


212±18 


219±88 


299±11 


< 150 


292±29 


223±132 


375±45 


< 189 


2 


117±15 


334±44 


226±13 


< 150 


213±21 


< 189 


203±36 


213±169 


1 


-4±9 


373±36 


116±8 


213±41 


152±18 


275±35 


160±27 


259±92 





-86±10 


496±34 


-4±10 


338±43 


32±18 


359±36 


62±33 


391 ±116 


■ 1 


-68±8 


383±25 


-49±9 


186±51 


-52±17 


273±34 


48±31 


304±111 


2 


-114±11 


171±54 


-84±9 


176±53 


-14±20 


238±52 


-95±31 


< 189 


3 


-156 ±10 


197±31 


-126±15 


203±82 


-112±19 


205±61 


-69±35 


313±136 


4 


-207±20 


347±87 






-5±42 


546±226 







a Positive position is southeast along the 160° slit. 

b Ccntroids arc relative to 2342 km s — 1 , the systemic velocity of NGC 2110. 
c The resolution has been subtracted in quadrature. 



Table 8: NGC 4388 Line Centroids and Widths 







[Fc II] 


Pa/3 




H 2 




Br-/ 




[S IXJ 


Pos a 
II 


Ccntroid fa FWHM C 
km s — ^ km s — ^ 


Centroid 
km s — 1 


FWHM 
km s — ^ 


Centroid 
km s — 1 


FWHM 
km s — 1 


Centroid 
km s — 1 


FWHM 
km s — ^ 


Centroid FWHM 
km s — 1 km s — 1 



NGC 4388, 30° slit 



4 


18±20 


353±60 


60±20 


192±75 






109±53 


454±237 






3 


35±37 


< 150 


45±16 


334±64 






35±33 


424±139 






2 


38±23 


443±94 


2±10 


206±30 


80±21 


197±114 


1±18 


284±71 






1 


-2±13 


249±42 


4±11 


168±40 


48±14 


338±53 


7±14 


286±52 


-70±17 


< 150 





-42±17 


479±67 


-12±12 


249±43 


8±10 


306±34 


6±15 


< 182 


-5±28 


< 150 


-1 


-51±12 


606±34 


-21±11 


366±48 


-60±11 


316=10 


23±35 


403±84 


-110±22 


204±95 


-2 


52±12 


624±31 


12±28 


321±62 


-63±19 


272±86 


-38±15 


392±53 






-3 


76±17 


195±64 


16±18 


383±70 


-50±27 


278±173 


158±29 


524±98 






-4 


138±22 


357±91 


24±27 


< 150 






















NGC 


4388, 120° slit 










2 


34±23 


601±103 


49±19 


218±81 


3±16 


253±80 










1 


10±21 


624±75 


91±20 


291±54 


2±10 


221±34 


56±17 


404±41 









28±13 


380±60 


54±5 


267±56 


-4±11 


211 = 10 


54±16 


286±62 


-29±11 


201±36 


-1 


14±17 


332±100 


47±6 


333±20 


-37±9 


284±30 


71±15 


331±58 






-2 


13±21 


252±140 


1±13 


240±67 


-42±12 


249±442 


67±34 


< 194 







a Positive position is towards the position angle of the slit as quoted. 

b Ccntroids arc relative to 2525 km s~ ^, the systemic velocity of NGC 4388. 

c The resolution has been subtracted in quadrature. 
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Table 9: Mk 3 Line Centroids and Widths 





[Fe 


II] 


P 


a/3 




H 2 


B 


r-y 


Pos a 


Centroid 6 


FWHM C 


Centroid 


FWHM 


Centroid 


FWHM 


Centroid 


FWHM 




km s" 1 


kms" 1 


km s" 1 


kms" 1 


km s" 1 


km," 1 


km^ 1 


km^ 1 


3 


19±26 


382±105 


73±12 


200±65 


72±60 


< 191 


24±28 


265±145 


2 


-44±16 


367±51 


-58±36 


255±125 






74±23 


< 191 


1 


-363±19 


1165±102 


-306±14 


1201±122 


-22±79 


710±335 


-153±46 


899±210 





-190±13 


1137±41 


-109±H 


991±46 


-2±49 


639±99 


-15±21 


564±169 


-1 


54±13 


501 = 56 


-4±13 


695±53 


81±34 


< 191 


127±19 


446±85 


-2 


162±16 


394±83 


-84±61 


532±220 


165±32 


745±430 


3±18 


< 191 


-3 


-77±45 


1117±248 


-75±21 


< 155 











a Positive position is east towards the 113° position angle of the slit. 
b Centroids are relative to 4046 km s~^. the systemic velocity of Mk 3. 
c The resolution has been subtracted in quadrature. 
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